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TECHNICAL  DISCUSSION 


I.  INTRODUCTION 

SAI  has  continued  to  achieve  notable  gains  in  under¬ 
standing  the  interactive  dynamics  of  radiative  emission, 
aesorption,  and  exchange,  and  plasma  hydrodynamic  motion. 
During  the  performance  period  of  the  contract  addressed 
by  this  report  (2"  August  1980  to  31  December  1981)  the 
r-dependent  collapsing  z-pinch  code  U'HYRAD  was  vastly 
improved  in  numerical  efficiency  and  flexibility  and  bench- 
marked  to  establish  its  accuracy.  The  code  was  applied  to 
various  experimental  conditions  and  elements  in  a  predictive 
and  analytical  effort  which  yielded  much  insight  into  the 
general  physics  as  well  as  specifically  valuable  design 
guidelines.  Important  advances  were  made  in  the  modeling 
of  the  radiation  itself  (apart  from  the  hydrodynamics), 
such  as  extension  of  escape  probability  techniques  to 

ft 

cylindrical  and  spherical  geometry.  Our  efforts  in  detailed 
radiation  modeling  are  covered  in  Section  2  and  Appendices  A 
through  C;  the  2-pinch  radiation-hydrodynamics  calculations 
are  described  in  Section  3  with  a  discussion  of  the  KHYRAD 
code  subroutines  included  in  Appendix  D. 


DETAILED  RADIATION  MODELING 


As  we  have  demonstrated  in  a  previous  NRL  contract 
and  in  publications,*'11  the  application  of  escape  probabilit 
concepts  to  line  transport  can  be  very  fruitful  in  reducing 
computational  labor  while  maintaining  an  adequate  degree  of 
accuracy.  Until  this  year  such  applications  have  been 
limited  to  plane  parallel  geometry  due  to  the  difficulties 
of  solid-angle  integrations  in  spherical  and  cylindrical 
geometry.  However,  we  have  discovered  a  method  to  bypass 
the  integrations  by  use  of  a  single  angle  for  inner-to- 
outer  cell  couplings  and  a  reciprocity  relationship  for 
coupling  constants  for  oute r- to  -  inner  cell  couplings. 

This  allows  the  extension  of  probabilistic  techniques  to 
cylindrical  and  spherical  geometry.  The  mathematics  is 
detailed  in  Appendix  A,  "Direct  Solution  of  the  Equation 
of  Transfer  Using  frequency-  and  Angle -Averaged  Photon- 
Escape  Probabilities  for  Spherical  and  Cylindrical 
Geometries . " 

Our  extensive  spherical  plasma  spectrum  simulation 
capability  utilized  for  NRL  is  aimed  at  employment  in 
pellet  experiments  where  the  actual  emitted  spectrum  is 
measured.  Before  such  an  approach  is  employed  with 

confidence,  the  usefulness  of  the  simulated  spectra  must 

(  ) 

be  benchmarked.  Last  year’s  report  presented  a  success¬ 
ful  one-on-one  comparison  of  a  calculated  neon  spectrum 
with  a  published  spectrum  obtained  from  a  pellet  implosion 


experiment  at  the  University  of  Rochester.  This  demonstrate d 
that  an  exact  replication  of  an  experimental  spectrum  from 
an  opticallv  thick  plasma  was  feasible  and  could  lead  to 
accurate  temperature  and  densitv  diagnostics. 

Since  then  a  great  many  more  pellet  spectrum  calcula¬ 
tions  have  been  performed  for  predictive  and  analytical 
purposes.  The  goal  of  this  latterset  of  calculations  is 
to  understand  how  the  spectrum  changes  with  pellet  tempera¬ 
ture  and  density.  The  calculations  are  detailed  in 
Appendix  B,  "Col  1 isional - radiat ive-equi 1 ibr ium  spectroscopic 
diagnosis  of  a  compressed,  optically  thick  neon  plasma." 

In  addition  to  elaborating  in  detail  the  expected  spectral 
changes  with  degree  of  heating  or  compression,  an  approximate 
analytic  model  of  the  approach  to  LTE  of  a  multistate, 
multistage  plasma  is  developed  in  this  appendix. 

Prediction  and  analysis  of  plasma  mixture  experiments 
and  the  effects  of  overlapping  absorption  and  emission  have 
been  of  continuing  interest  to  SAI  and  NRL .  Ke  have 
recently  considered  photopumping  in  opacity  resonances  for 
Silicon/Aluminum  and  Sodium/Neon  mixtures.  Our  detailed 
results  are  given  in  Appendix  C,  "Plasma  Conditions  Required 
for  Attainment  of  Maximum  Gain  in  Resonantly  Photo-Pumped 
A 1  XII  and  Xe IX." 

Until  this  past  year,  our  analyses  of  wire  plasma 
spectra  have  concentrated  on  inferring  average 


temperatures  and  densities  for  these  plasmas.  To  this  end, 
spectra  from  assumed  isothermal,  isodense  plasmas  have 
been  computed  and  compared  with  the  measurements  in  an 
iterative  process  until  agreement  was  obtained.  This  was 
satisfactory  as  an  initial  procedure  in  obtaining  approxi¬ 
mate  characteristic  plasma  temperatures  and  densities. 
Moreover,  an  understanding  of  more  realistic  situations 
where  gradients  are  present  may  only  be  obtained  if  the 
simpler  homogeneous  plasmas  and  the  evolution  of  their 
spectra  are  fully  comprehended.  This  year  we  have  extended 
our  previous  first  order  analyses  of  homogeneous  plasmas 
to  investigate  the  effects  of  gradients. 

Some  of  the  ambiguities*'11  which  may  result  from 
using  line  ratios  only  for  spectroscopic  diagnostics  may 
be  at  least  partially  resolved  by  use  of  hi gh-d i spers i on , 
high  resolution  spectroscopy  in  analyzing  line  prof iles 
along  with  their  measured  intensity.  Temperature  gradients, 
in  particular,  are  susceptible  to  diagnosis  using  line 
profiles,  since  the  degree  of  sel f- reversal  in  an  optically 
thick  line  depends  on  the  radial  temperature  profiles.*"*1 
A  full  physical  explanation  for  this  effect  appears  in 
Reference  4.  As  detailed  in  our  final  report  for  last 
year,*"1  plasma  implosion  and  its  accompanying  frequency 
Peppier  shifts  may  transform  a  symmetric  line  profile  self- 
reversal  into  an  asymmetric  one,  with  the  blue  wing 
dominant.  In  principle,  therefore,  both  implosion  velocity 


4 


and  temperature  profile?  may  be  inferred  by  analyzing  an 
optically  thick  profile.  he  present  the  following  analyse? 
of  exploded- wire  spectra  as  successful  applications  of  the 
usefulness  of  these  effects. 

A  first-order  analysis  of  the  spectrum  of  PI  shot  329 

has  been  presented  elsewhere . 1 ^ !  In  the  present  analysis 

we  utilize  the  approximate  information  obtained  previously 

to  infer  additional  plasma  properties  from  the  Si  XIII 
]_ 

Is- -Is  2p  I'  resonance  line  profile.  As  pointed  out  in 
reference  5,  this  line  seemed  to  exhibit  opacity  broadening, 
with  a  peculiar  "shoulder"  present  on  the  red  wing  of  the 
line.  A  high  dispersion  spectrum  of  the  appropriate  wave¬ 
length  region  appears  in  fig.  1.  If  the  average  conditions 

n 

diagnosed  in  ref.  5  (.1  =  050  cV,  n  =  2  x  10“  .)  are  applied 
in  modeling  the  resonance  line  profile  with  an  implosion 
velocity  of  4  x  10(1  cm  sec  1  at  the  outer  plasma  surface, 
one  obtains  from  a  mu  1 1 i f requency  calculation  the  resonance 
line  profile  plotted  in  fig.  2.  This  profile  was  computed 
assuming  an  isodense,  isothermal  plasma  and  a  pinhole- 
measured  radius  of  250  um  (P.  Burkhalter,  private  communica¬ 
tion;.  Note  the  striking  resemblance  to  the  measured 
spectrum  of  fig.  1.  However,  once  the  resonance  line  pro¬ 
file  is  convolved  with  the  estimated  3  m*  of  experimental 
broadening  (fig.  3)  the  shoulder  feature  washes  out  and 
the  resemblance  disappears.  The  isodense,  isothermal 
model  therefore  does  not  work  when  the  finite  resolution 
of  the  spectrograph  is  properly  accounted  for. 


I  ho  results  o!  a  vast!  v  improved  profile  s  iir.u  lat  i  c-r. 
for  this  shot  are  presented  in  figs.  4  and  5.  figure  - 
presents  the  computed  line  profile  ■ wi theut  experimental 
broadening  added  in  for  a  plasma  identical  to  that  of 
fi.es.  2  and  5  e .\ c e p t  that  the  temperature  profile  decreases 
sharply  as  radius  increases.  he  inner  201  of  the  plasma 
is  assumed  to  have  a  temperature  of  11“0  eV ;  from  20‘,  to 
1001  of  the  total  radius  'I  is  presumed  to  ramp  linearly 
downward  from  500  to  1 ' 0  eY.  This  cold  plasma  exterior 
serves  tc  improve  the  profile  agreement  by  producing  a 
deep  sel f- reversal  (sec  fig.  4i.  Even  after  convolution 
with  a  5  mA  instrumental  gaussian,  the  profile  is  in 
excellent  agreement  i  fig.  5>  with  experiment.  This, 
inference  of  a  temperature  gradient  is  fully  in  accord 
with  our  new  analysis  of  another  shot  which  we  now  discuss. 

Based  on  an  investigation  of  the  effects  of  gradients 
on  the  emitted  spectra  of  plasmas  by  some  of  the  NRL  staff' 
we  have  broadened  and  expanded  our  analysis  of  Maxwell 
Laboratories  shot  159.  This  analysis  further  supports  the 
inference  of  temperature  gradients  of  the  type  described 
above  in  PI  shot  529. 

Shot  159  resulted  from  the  implosion  of  A1  wires  and 

the  measured  t i me  - i n tegrated  spectrum  exhibits  an  asymmetri 

•>  ] 

sel f- reversed  prolile  lor  the  A1  XII  Is-- Is  2p  P  resonance 
line  | fig.  6).  Symmetric  sel f- reversed  profiles  are  a 


0 


wel 1  -  unde rs t ood ,  classical  effect  of  large  line  enacitv. 
lo  produce  asymmetry  with  a  suppressed  red  u  my,  however, 
the  Poppl er  shift  of  plasma  implosion  combined  with  large- 
line  opacity  is  required,  as  detailed  in  last  year's  report, 
he  previously  obtained  an  excellent  fit  to  this  profile- 
assuming  a  homogeneous  plasma  of  radius  250  ,  temperature 

of  500  eV ,  ion  density  s  x  10^  cm  ' ,  and  implosion  velocity 
4  x  1 0°  cm  sec  l.  However,  this  fit  required  an  experimental 
spectral  resolution  s  1 .  t>  m  in  order  that  the  profile 

not  wash  into  a  gaussian  due  to  experimental  broadening. 

This  resolution  is  a  factor  of  2-5  better  than  the  estimates 
(.which  are  themselves  somewhat  uncertain!.  The  washout  to 
a  gaussian  could  be  prevented  if  the  intrinsic  profile  (the 
one  which  would  be  measured  with  infinite  resolution!  were 
wider  with  a  deeper  core  scl f- reversal .  Such  a  profile 
can  be  produced  by  an  i nhomogeneous  plasma  with  a  hot  core 
and  colder  outer  region.  he  have  found  that  by  employing 
such  a  gradient  in  our  simulation  that  the  measured  profile 
can  be  reconciled  with  larger  experimental  broadening  and 
thus  be  completely  consistent  with  all  available  information. 
A  slightly  smaller  velocity  of  2  x  10*'1  cm  sec  1  is  also 
required.  This  velocity  is  interpreted  as  the  mean  velocity 
of  implosion  during  the  K-sholl  emission  pulse.  he  have 
also  simulated  the  t ime- i ntegrat i on  of  the  spectrum  by 
summing  the  line  profiles  produced  by  a  range  of  implosion 
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velocities  from  0  to  4  x  10L’  cm  sec  '  ,  assuming  the  velocity 
changes  linearly  with  time.  The  result  is  that  the  resultant 
"integrated"  theoretical  profile  is  virtually  identical  to 
that  produced  hy  snapshot  emission  from  a  plasma  imploding 
at  I  x  10h  cm  sec  Therefore,  the  velocity  ref lected  by  the 

shape  o  f  the  em  i  t  ted  1  i  ne  profile  is  consistent  with  the  physical 
mean  velocity  of  implosion.  This  is  to  be  expected,  as  a 
zero  velocity  plasma  will  produce  a  symmetric  profile  and 
one  of  say  4  x  1 0*’  cm  sec  ^  will  produce  a  profile  more 
skewed  than  at  1  x  10°  cm  sec  *.  1'he  sum  of  the  emitted 
profiles  will  thus  be  characteristic  of  an  intermediate 
velocity  which  is  also  consistent  with  a  physical  mean. 

Of  course,  to  produce  this  profile  there  is  no  escaping 
the  conclusion  that  the  bulk  of  the  K-shell  emission  must 
occur  during  implosion  rather  than  expansion  after  the 
peak  compression. 

The  precise  conditions  we  have  inferred  from  line 
ratios  as  well  as  the  resonance  line  profile  are  as  follows: 
the  inner  20','  •  (radiallv)  is  at  900  eV ;  the  outer  80i  is 
at  400  eV  just  outside  the  900  eV  core  and  ramps  linearly 
to  1x0  eV  at  the  outer  boundary  of  the  emitting  plasma. 

As  is  shown  in  the  z-pinch  solutions  given  in  Sec.  3  the 
region  of  appreciable  radiative  yield  is  strongly  depend¬ 
ent  on  the  plasma  density.  In  low  emission  regions  with 
the  high  ionization  level  obtained,  corona  temperature 


S 


can  rise  to  It'  koV  or  more,  but  direct  confirmation  is 
lacking.  The  total  ion  density  remains  unchanged  from  our 
previous  analysis  at  8  x  ll1^  cm  The  plasma  diameter 

of  1  mm  is  roughly  consistent  with  the  somewhat  irregular 
configuration  revealed  by  the  pinhole  picture.  The  profile 
produced  by  such  a  plasma  is  presented  in  figs.  ~  through  9 
for  assumed  instrumental  linewidths  of  0,  2,  and  3  m.k , 
respectively.  Note  that  the  experimental  smearing  of  the 
intrinsic  profile  of  fig.  “  gradually  makes  the  self¬ 
reversal  shallower,  to  t  lie  point  where  fig.  9’s  profile 
nicely  matches  experiment.  The  assumed  resolution  of  3  nu 
is  in  excellent  agreement  with  estimates  based  on  the  exper 
mental  con  f  i  gura  t  i  on  .  figure  11'  presents  the  profidc  of 
the  resonance  line  which  would  be  produced  by  an  identical 
s  t  a  t  i  c  plasma  (no  implosion).  This  line  profile  bears  no 
resemblance  to  the  experimental  onc--which  fact ,  combined 
with  the  good  resemblance  of  fig.  9--is  clear  evidence  for 
implosion  effects  on  the  line  emission  in  the  laboratory 
frame.  Ihe  character  of  the  inferred  temperature  gradient 
closely  resembles  that  of  the  silicon  shot  329.  Taken 
together,  these  two  analyzed  spectra  demonstrate  the  power 
of  spec t roscop i e  analysis  in  inferring  most,  if  not  all, 
of  tlie  plasma  properties  of  interest. 
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FIGURE  CAPTIONS 


Fig,  I 

Fig.  : 

Fig.  3 

Fig.  4 

Fig.  5 

Fig.  b 

Fig. 


">  1 

Fligh  dispersion  spectrum  of  the  SI  XIII  ls“-ls  2 p ^ P 
and  intercombination  lines  obtained  from  PI  shot  329. 
Theoretical  resonance  and  intercombination  line 
profile  obtained  from  a  250  um  Si  plasma  of 
implosion  velocity  4  x  10b  cm  sec  Assumed 

•’0  .  i 

temperature  is  050  eV;  electron  density,  2  x  10“  cm  '. 
Same  as  Fig.  2,  except  that  the  spectrum  has  been 
convolved  with  a  5  m  „  gaussian  to  simulate  finite 
experimental  resolution. 

Silicon  plasma  spectrum  calculated  to  arise  from 
the  same  plasma  modeled  in  Fig.  2,  except  that  a 
temperature  decrement  (described  in  the  text  i  has 
been  applied. 

Same  as  Fig.  4,  except  that  the  spectrum  h3s  been 
convolved  with  a  5  m  gaussian  to  simulate  finite 
experimental  resolution. 

High  dispersion  spectrum  of  the  A1  XII  principal 
resonance  and  intercombination  lines  obtained 
from  Maxwell  Laboratories'  shot  159. 

Computed  spectrum  which  would  be  emitted  by  a  1  mm 
diameter  A1  plasma  of  ion  density  8  x  10^  cm 
The  temperature  profile  is  as  described  in  the 


1 1 


text.  The  implosion  velocity  is  assumed  to  be 
2  x  10b  cm  sec  \  .No  experimental  gaussian 
convolution  was  performed. 

Fig.  8  Same  as  Fig.  ",  except  2  m  A  experimental 
broadening  is  assumed. 

Fig.  9  Same  as  Fig.  ",  except  5mA  experimental 
broadening  is  assumed. 

Fig.  10  Same  as  Fig.  9,  except  a  zero  implosion  velocity 
for  the  plasma  is  assumed. 
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RAP  I  AT  I  OX  HYDRODYNAMICS-  COMI’AKA'I  I  YE  TRANS  I LNI 
CALCULATIONS  OF  IMPLODING  LOADS 


Section  1 
INTRODUCTION 

It  is  generally  accepted  that  cylindrical  geometry 
is  the  most  realistic  simple  geometry  for  modeling  col¬ 
lapsing  c-pinches.  Substantial  effort  has  been  expended 
on  the  establishment  of  a  1-D  radial  collapsing  c-pinch 
hydrodynamics  code  1YHYRAD,  which  can  model  experiments 
initialized  by  wires,  gas  puffs  or  annular  foils,  employing 
a  variety  of  ionization  physics  and  radiation  subroutines. 
This  code  is  a  significant  improvement  of  an  earlier  model, 
WHY RAC .  ^  A  variety  of  initial  radial  distributions  of 
particle  number  density  and  temperature,  temporal  variation 
of  machine  discharge  voltage,  and  anomalous  resistivity 
functions  may  be  employed.  Operational  flexibility  due  to 
general  re-start  and  improvement  of  computational  vectcriza- 
tion  was  also  incorporated.  This  report  summarizes  the 
analytical  and  numerical  assumptions  contained  in  the  code 
development,  illustrates  subset  benchmarks  made  to  prove 
modeling  adequacy,  and  applies  the  model  to  various  aluminum 
and  argon  experimental  conditions.  Radiation  yields  and 
power  are  given  in  normalized  magnitudes  so  as  to  maintain  an 
unclassified  status;  direct  yield  values  are  given  in  a 
companion  DNA  final  report  which  does  not  include  some  of 
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the  benchmark  and  test  case  deve 1 opment .  Appendix  D  is 
also  included- in  this  unclassified  version  which,  documents 
the  code  structure  by  describing  the  function  of  the  variou 


subroutines . 
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» 


» 


♦ 


t 


> 


i 


* 


PLASMA  DYNAMIC  MODEL 


Analysis  of  the  dynamics  of  a  radiation  coupled 
imploding  :-pincn  plasma  depends  on  solutions  of  the  time- 
dependent  equations  of  mass,  momentum,  electron  tempera¬ 
ture  and/or  total  internal  energy,  magnetic  flux  and  the 
equations  of  state,  as  well  as  the  specification  of  appropri¬ 
ate  transport  processes.  The  solution  of  the  equations 
are  linked,  at  every  step  and  radial  cell,  to  a  radiative 
flux  model  in  which  both  line  and  continuum  emission  and 
opacity  are  controlled  by  a  probability  of  photon  escape 
representation  of  cylindrical  radiation  transport  solutions. 
Solution  of  these  equations  is  also  dependent  on  the  coupling 
with  the  driving  discharge  circuit  which  in  turn  is  altered 
by  the  changing  plasma  induction.  The  energy  density  of 
the  plasma  involves  an  adiabatic  exponent  containing  the 
potential  ionication  energy;  charge  neutrality  is  also  assumed. 

The  WHY RAD  System  solves  the  following  fluid 
conservation  equations: 
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Momentum : 
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Electron  Temperature: 
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( 4 ) 


Ion  Temperature  for  Total  Internal  Energy  Balance 
(optional j  : 


(5) 


The  dependent  variables  in  the  above  conservation  equa¬ 
tions  are  p,  plasma  mixture  density  (=>-  ion  density);  y, 
radial  plasma  velocity;  p,  plasma  pressure;  i,  axial  cur¬ 
rent  (c-direction) ;  B,  atimuthal  (e -direction)  magnetic 
field;  E,  total  plasma  internal  energy;  e,  axial  electric 
field;  total  radiative  power;  ne ,  p  ,  Tp ,  electron 

number  density,  pressure,  and  temperature,  respectively. 
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Qe  and  qg  are  non- radiative  thermal  flux  terms,  given 
below.  The  adiabatic  exponent,  ■>  ,  which  is  a  function  of 
electron  temperature,  number  density  and  ionization  number 
{see  Eq.  (5")},  is  assumed  uniform  ever  an  individual 
spatial  cell  or  time  step  such  that  its  gradients  are 
neglected.  It  is  conceivably  possible  that  a  stable 
iterative  solution,  including  the  Tg  ,  ng>  and  Zj  gradients 
in  the  y  expression,  within  each  cell  and  time  step,  may 
be  constructed:  however,  D.  Colombant  (private  communica¬ 
tion)  found  in  work  with  the  h'HYRAC  code  that  numerical 
instabilities  seemed  to  be  quickly  generated  with  this 
modeling.  In  addition,  computational  run  times  would 
rapidly  increase  so  that  practical  considerations,  forced 
the  abandonment  of  this  refinement. 

The  Qe  term  in  Eq.  (4)  is  given  by 


am. 


m. 
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(6) 


where  is  the  electron-ion  equilibration  time.  For  the 

general  non- tempe rature  equilibrated  situation  incorpora¬ 
tion  of  the  second  term  in  Eq.  (6)  is  performed  by  a 
subsequent  Te  evaluation  at  the  end  of  each  time-step 
cycle,  i.e.,  a  solution  of 
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This  relationship,  Eq.  ("),  represents  a  common  scheme 
employed  in  explicit  convective  hydrodynamic  codes;  cer¬ 
tain  flux  terms  have  to  be  treated  in  a  time- split  manner 
and  solved  after  the  convective  solution  for  T  ,  so  that 
one  updating  of  Te  within  a  time  step  involves  the 
separate  solution  of  Eq.  (7).  A  solution  of  Eq .  (7)  for 

T  in  terms  of  the  exnlicitlv- obtained  T  and  the  ion 
e  e 

temoerature  Tt  is  given  in  Ref.  1:  however,  T.  values 
computed  from  Eq.  (5)  seem  to  be  excessively  high  indicat¬ 
ing  that  future  calculations  should  employ  a  full  twc- 
temperature  pair  of  energy  equations  including  the 
equilibration  term,  Eq .  (6).  Almost  all  the  solutions 

illustrated  in  this  report  assume  equal  electron  and  ion 
temperatures  so  that 


The  q£  electron  thermal  flux  term  is  (ion  thermal  fluxes 
are  neglected)  given  by: 


where 


and  £  are  given 
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is  given  by  Eq .  (40a) ,  which  is  employed  equivalently 
to  in  Ref.  2.  The  6,  y',  S"  coefficients  are  func- 

f  *>  i 

tions  of  Zj,  the  average  plasma  ion  number.  Braginskii 
gives  their  magnitudes  at  five  2j  values  between  1  and  « 
quadratic  fits  have  been  constructed  which  are  within  a 
few  percent  accuracy  for  Zj  >  5.  The  results  are: 


5  =  0.31652 
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-J 
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S'  =  1.282 

-  0.0242 

2,  -  5.600 

x  10*4 

"I 

S'  =  1.5 

•>  '  =  5.350 

-  0.17855 

2  j  +  5.66" 

x  10' 

“  I 

y':  =  3.622 

-  0.0222 

2j  +  3.600 

x  10"  4 

"1 

The  axial  current  velocity  (drift 

steed) 

is 

given  by 

U_  =  -j/eng  .  (10) 
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Thermally  perfect  state  expressions  for  electron  pressure 
and  total  pressure  and  the  charge  neutral  plasma  electron 
density  are  given  by 
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Equations  (Ij  through  (4,,  not  including  the 

r  ^SfrKp  ST  /Sri 'Sr  thermal  diffusion  and  the  temperature 
^  x  e 

equilibration  terms,  are  solved  by  the  explicit  Flux- 

f  ~  1 

Corrected- Transport  algorithm  ETBFCT.  ' 1  The  plasma  is 
allowed  to  flow  between  the  individual  cells,  but  the 
outer  boundary  of  the  grid  conforms  to  the  displacement 
specified  by  a  mean  velocity  over  the  high  density  cells, 
ever  that  time  step.  The  details  of  the  outer  boundary 
displacement  bnd  other  aspects  of  the  reconing  algorithm 
are  given  at  the  end  of  this  section.  The  FCT  method 
solves  a  general  continuity  relation  of  the  form 
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where  the  last  three  terms  on  the  right  hand  side  would 
model  a  variety  of  source  function  representations .  The 


algorithm  is  basically  an  imp rovement  of  the  Lax-Kendrcf f 
technique  whereby  strong  shocks  or  high  fluid  gradients 
are  preserved  and  appropriately  convected  through  the 
fluid.  It  adds  a  diffusion  term  during  convective  trans¬ 
port,  and  subtracts  it  out  over  most  of  the  flow  field  in 
the  antidiffusive  phase  (second  half!  of  each  time  step. 

Any  residual  diffusion  is  sufficient  to  prevent  dispersive 
ripples,  and  ensure  positivity  of  the  fluid  density  func¬ 
tion  t .  The  ETBFCT  algorithm  is  not  restricted  to  a  strict 
fixed  grid  Euierian  setup  or  uniform  cell  widths,  so  no 
basic  coding  difficulty  is  present,  although  too  large  a 
change  of  the  mesh  over  a  time  step  will  destroy  the 
accuracy  of  the  algorithm.  Equations  (1)  through  (4)  have 
convective  densities  of  c ,  : v ,  E,  and  T0  while  the  source 
terms  describing  pressure  gradient,  Lorentr  force,  Joule 
heating,  radiation  transport,  and  magnetic  field  dependent 
thermal  fluxes  fit  into  the  above  forms. 

Incorporation  of  the  Ke^?Te/3r  part  of  the  thermal 
conductivity  is  also  by  a  time-split  evaluation  at  the  end 
of  each  time-step  cycle,  i.e.,  after  successive  calcula¬ 
tions  of  Tg  employing  convective  FCT  and  electron-ion 
equilibration  (when  desired),  then  the  thermal  diffusion 
influenced  Tg  is -obtained  by  fully  implicit  solution  of 
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using  a  tridiaeonal  method.  Here  fT  is  the  addition  tc 
-  e 

the  FCT  and  equilibration  results.  Comparison  of  Eqs.  ',5 
and  (41  show  that  the  related  internal  energy  change  is 
given  as 


£E  =  |r.e  (•>-!)}  f(k  Te/ 

(with  Te  in  eY  and  E  in  erg/cnTl.  As  indicated  previously, 
the  adiabatic  exponent  is  assumed  uniform  over  a  cell  and 
time  step  such  that  gradients  of  its  functional  dependency 
on  Tg ,  ngand  I T  are  neglected.  This  assumption  facilitates 
the  above  proportionality  between  ;E  and  fTe  although  it 
is  an  approximation. 

Before  considering  the  boundary  conditions  on 
these  conservation  equations  it  is  necessary  to  specify 
the  magnetic  flux  induction  equation.  The  3B/ct  =  -cv  x 
Faraday  equation  is  combined  with  a  generalized  Ohm’s  law 
relating  the  electric  field  to  electric  current,  Lorentz 
force,  pressure  and  temperature  gradient  contributions. 

It  is  assumed  that  the  induced  field  is  represented  by 
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terns  incorporate  magnetic 


thermal  diffusion  contribution.  The  function  r.  is  an 

assumed  resistivity.  Solution  of  this  relation  is  also 

by  a  fully  implicit  tridiagonal  algorithm  cue  to  E-field 

propagation  swiftness  with  non- zero  values  of  r, .  The 

T  distribution  used  in  Eo .  ,151  is  the  result  from  the 

e  ■ 

previous  time  step  as  the  B  field  has  to  be  computed 

prior  to  the  hydrodynamic  conservation  equations;  as  no 

back  and  forth  iteration  solution  between  B  and  T  was 

e 

attempted  it  was  deemed  necessary  to  incorporate  a  more 
stringent  time  step  criteria  employing  a  restrictive 
change  to  the  total  internal  energy,  Bq .  (25d).  Iterative 

close  coupling  between  Tg  and  B,  to  ensure  that  they 
propagate  exactly  together  is  desirable  but  the  computa¬ 
tional  time  penalty  for  large  cell  number  hydrodynamics 
was  not  considered  warranted  for  either  U'HYRAC  or  WHY  RAD . 

The  boundary  conditions  on  the  system  of  equations 
are  (R  is  the  outer  boundary  of  the  implosion): 
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Establishment  of  the  radiallv  intecrated  tctal  discharge 


current  l(t)  yields  the  outer  boundary  condition  on  the 
E  field,  leading  to  its  solution  and  the  specification  of 
the  Lcrento  force  anc  Joule  heating  source  functions. 
Therefore,  appropriate  coupling  into  the  external  driving 
circuit  is  crucial;  denoting  the  external  circuit  impe¬ 
dance  by  1^.  inductance  by  Lc,  driving  voltage  by  Y  (t) 
and  effective  voltage  across  the  plasma  by  Y  ft)  we  have  that 


Yp(tj  is  obtained  by  integration  cf 


T  X  £  =  - (l/c . (cB/ct; 


or 


o  £  _  1  c  B 

cT  c  Tt 


(16) 


(16a  j 


across  a  loop  along  a  perfectly  conducting  outer  shell 
(at  a  radius  r  >>  R(t)}  and  cathode  end  plate,  through 
the  plasma  at  a  general  radius  r,  and  across  the  Yp  voltage 
drop  at  one  end  of  the  discharge  (see  Figure  1>.  The  axial 
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where  Yp,'  the  plasma  voltage,  is 

YpU,i;  =  -  vC  (c 

o 

Here  r  is  the  outer  radial  location  of  the  perfectly 
\\ 

conducting  return  current.  The  function  V  (r,t)  can 
equivalently  be  written  as 
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where  the  local  electric  field  is 

T* 

e(r)  =  n(o)i(c)  +  |  J  f|  dr'  •  (20d 

o 

This  expression  for  e  is  consistent  with  a  radial  integra 
tion  of  the  Faraday  law,  Eq.  (16a),  together  with  use  of 
Eq.  (1")  at  r  =  0  and  the  B|r_Q  =  0  boundary  condition. 
Equation  (20c)  is  a  representation  of  the  plasma  voltage 
around  the  loop  in  Fig.  1.  The  major  electrodyanmic 
assumption  in  the  W’HYRAC^  model,  and  retained  with  some 
misgivings  in  the  WHY RAD  code,  is  that  a  single  value  of 
at  the  cell  location  where  the  current  stops  flowing  for 
the  previous  time  step  (r  =  R) ,  is  employed  to  calculate 
Yp .  Then  Eq .  (20c),  can  be  solved  to  give 
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I/Rc  so  that 
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e-R)  is  Eg.  (1“)  evaluated  at  r  =  R  and  at  the  previous 
time  step, 
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with  the  current  cutoff  at  r  =  R.  It  was  originally 
considered  in  Kef.  1  to  directly  employ  Eq.  ( 1 ” j  in  the 
calculation  of  t '  r 1  for  all  r  '.not  i us t  at  r  =  R ;  .  How¬ 
ever,  the  resulting  field  was  very  choppy  cue  to  the 
influence  of  large  variations  of  velocity,  number  density, 
and  temperature.  The  tost  way  to  calculate  c(r)  is  from 
direct  radial  integration  of  the  Faraday  law,  Eo.  (16a), 
after  specification  of  the  B-field  vEq.  20d.l  .  Insertion 
of  E q .  (21)  into  E q .  (15)  g i v e s 


Vt; 


LcClK 


whi  ch 


is  numerically  integrated  in  time  to  specify  the 


tc-tal  current  and  the  outer  boundary  condition  on  the 
E  field. 

The  local  current  is  obtained  from  the  Maxwell  law 
in  the  MHD  approximation 


which  is  solved  by  a  central  finite  difference  representa¬ 
tion  of  the  gradient;  as  indicated  above  the  e  field 
distribution  is  then  obtained  by  employment  of  Eq .  (20d). 
The  electromagnetic  field  dependent  source  terms  in  the 
hydrodynamic  equations  are  therefore  determined. 

Time  step  criteria  is  evaluated  from  the  cell- 
dependent  minimum  of  the  following  relations. 
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Joule  Heating  Energy  Transfer: 
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Radiative  Energy  transfer : 


ct  =  min 
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Previous  Time  Step  Internal  Energy  Change : 
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Direct  incorporation  of  the  Alfven  speed,  =  ^/b"/4*:c  , 

into  a  time  step  criterion  such  as  Eq .  (25a)  is  unwise  as 

the  presence  of  a  large  B  field  together  with  a  low  density 

in  the  outer  regions  of  the  implosion  yields  excessively 

small  time  steps.  In  this  situation  the  time  steps  are 

chosen  at  a  location  which  is  far  removed  from  the  high 

density,  high  radiative  regions  of  interest.  The  WHY RAC 
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scheme  of  Ref.  1  is  retained,  wherebv  a  term  B"/4r(Cti“ 

A ' 

is  added  to  the  inertia  in  the  low  density  region.  is 

an  effective  Alfven  speed  of  reduced  magnitude  determined 
from  the  assumptions  that 
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t  from  the  other  time  step  criteria, 
d  t  . 
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This  ensures  that  the  Alfven  time  scale  is  no  faster 
than  the  hydrodynamic  time  scale.  Even  though  this  effect 
of  the  Alfven  speed  is  included  in  the  formalism  it  was 
decide'd  in  the  construction  of  WHY  RAD  to  include  the 
internal  energy  change  criteria,  Eq.  ( 2 5 d ) ,  in  addition 
to  the  first  three  criteria,  thus  practically  yielding  a 
significantly  finer  time  step  resolution.  This  last 
criteria  was  not  included  in  h'HYRAC.  The  absolute  minimum 
over  all  the  cells  of  the  four  functions,  Eqs.  (25a) 
through  (25d;  is  selected,  with  time  steps  of  20%  to  50% 
of  the  minimum  Courant  value  normally  computed. 

For  no  retoning  the  spatial  increment  widths  are 
an  assigned  fraction  of  the  total  beam  radius  R,  defined 
by  ( R  =  RADR) 

RADR  *  RADRO  +  YGM  x  ft 1/2  (26) 

where  RADRO  is  the  old,  previous  time  step  value  of  the 
1  / 

total  radius,  ft  is  the  half  time  step  (ETBFCT  is 
applied  twice  successively  ever  a  full  time  step)  and  YGM 
is  a  velocity  given  by 

A,.  (26a) 
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where  i  is  the  cell  number,  X^-  is  the  number  of  cells  in 
the  summation,  and  v.  is  the  individual  cell  velocity.  The 
summation  is  not  carried  out  over  all  cells,  but  only  ever 
those  cells  with  a  density  at  least  IDv  of  the  maximum 
density,  which  gives  a  weighting  to  the  cells  with  maximum 
mass  flux  in  the  evaluation  of  the  beam  contraction. 

A  density-dependent  reconing  option  is  also  present 
in  the  code  which  assigns  the  minimum  cell  width  at  the 
maximum  density  and  vica  versa,  thus  increasing  the  cell 
density  at  regions  of  maximum  radiation  yield.  Kith  no 
re  coning 

Lri  =  RADR/XR  =  R/XR  (26b) 

where  \R  is  the  total  number  of  cells.  Although  more 
detailed  consideration  of  the  required  cel  1  number  is  given 
later  in  this  report  it  should  be  noted  that  around  40  should 
be  the  minimum  number  to  represent  the  i  current  and  E  field 
radial  distributions.  The  reconing  option  uses  Eq .  (26b) 

to  define  a  representat ive  uniform  cell, 


Ar  =  RADR/XR  , 


together  with  an  average  number  density  given  by 


( 2  ~a 


XR 

n°  ‘IE  ‘“i  '■  Vt.,t|/R 

i  =  1 


RADR  . 


(  2  "  b  _ 
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Then  the  retoning  is  dictated  by 


*  .0,0,  ,  M 

dr?  -  dr  (n  /n^l 


l::c} 


where  the  power  M,  to  be  careful,  should  be  <<  1.  Now. 


NR 


i  =  l 


RAD NEK 


which  is  not  necessarily  the  same  as  RADR ,  previously 
obtained.  Therefore,  the  dr-’s  are  rescaled  by 

ir.  =  dr?  x  (RADR/RADNEK;;  (27d) 

which  is  the  chosen  retoned  cell  width.  In  either  cell 
distribution  algorithm  the  cell-center  locations  are 
determined  from 

r.  =  0.5  dr ,  (28a) 

1  i 


and  r^  =  r ^ _ 1  *  0.5(dr.  *  dr^^)  »  i  >  1  .  (28b) 

I 

It  is  necessarv  to  constrain  the  M  tower  in 
Eq .  (27c)  to  prevent  retoning- induced  numerical  instabil¬ 
ities.  The  function  employed  is 


M 


M°  (1  -  exp  L - D] } 

f _ * -  . 


(29) 


where  D  *  min 


(-9a) 


the  b  anc  c  coefficients  have  been  determined,  by  numeri¬ 
cal  experimentation,  to  be  of  order  10  ",  If  n° |  refers 
to  the  absolute  value  of  the  change,  over  a  time  step,  of 
the  logarithm  of  the  average  number  density  given  by 
Eq.  (27b),  and  RE. is  a  function  of  the  local  energy 
balance  ratio.  Although  more  will  be  said  later,  an  M° 
value  of  0. 03-0.05  was  employed  in  test  cases.  The  energy 
balance  ratio  is  the  plasma  internal  energy  at  t  =  0  plus 
the  electromagnetic  energy  added  up  to  the  time  in  ques¬ 
tion  minus  magnetic  field,  emitted  radiation  and  boundary 
flux  energies,  all  divided  by  the  plasma  internal  energy 
at  that  time.  For  complete  energy  balance  this  ratio 
should  be  unity.  For  the  ratio  >  1 ,  RE  *  1/ratio-,  for 
the  ratio  <  1 ,  RE  =  ratio,  while  for  complete  energy 
balance  the  contribution  of  the  second  term  in  Eq .  C29a) 

is  neglected.  The  two  damping  contributions  found  in 
Eq.  (29a)  control  the  effects  of  the  sudden  number  density 
and  energy  balance  changes  on  the  retoning  algorithm. 
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Section  III 

RADIATION,  IONIZATION  AND  TRANSPORT  PROCESSES 

Three  separate  radiation- ioni cation  packages  have 
been  employed  in  the  implosion  solutions;  the  first  one  is 
ATFHYS,  a  relatively  simple  ground  plus  one  excited  state 
ion  balance  Kith  approximate  radiation  transport,  TERRYS, 
a  curve  fit  representation  of  the  icni cation- radiation 
package  MCELL  solutions  with  radiation  transport  assumed 
to  be  a  double  path  probability  of  escape  function  out  of 
the  plasma  volume  from  each  individual  radiating  cell,  and 
real  time  application  of  MCELL.  MCELL  is  an  ioni cation 
collisional- radiative  equilibrium  package  which  solves  for 
the  ground  and  excited  state  population  densities;  line, 
radiation  recombination  and  Bremsstrahlung  radiation  trans¬ 
port  are  modeled  with  emission  and  absorption  balanced  out 
ever  the  multi-cell  radiation  increment  cones. 

The  ATPHYS  modeling  of  argon  was  useful  to  estab¬ 
lish  the  scheme  of  running  the  code,  remove  the  numerous 
bugs  and  hardwired  quantities  present  in  the  original 
V.'HYRAC  source  deck,  and  incorporate  a  number  of  the  fea¬ 
tures  of  the  NHYRAP  development.  No  numerical  results 
from  ATPHYS  calculations  will  be  presented  in  this  report, 
however.  Real  time  employment  of  MCELL  was  the  original  i 
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intent  c£  the  c omputat i cns ;  however,  calculations  with  a 
many- line  spectra  aluminum  model  showed  that  in  some  cases 
9Si  of  the  run  time  per  time  step  was  taken  up  by  the 
radi at  ion- ioni cation  package  and  only  Z\  by  the  remaining 
hydro- induced  B  field-driving  circuit  algorithm.  Even 
though  this  setup  probably  represents  the  ideal  modeling 
of  the  nonequilibrium  coupled  hydrodynamic -radiation 
scheme,  it  was  felt  that  running  costs  as  well  as 
an  excessive  number  of  restart  inputs  over  a  typical 
3000  time  step  case  prevents  the  present  implementa¬ 
tion  of  the  full  real  time  scheme.  It  is  suggested  that  an 
energetically  optimum  MCELL  code,  tailored  to  run  10  times 
faster,  would  facilitate  this  highly  desirable  direct 
coupling.  Nevertheless,  post -processing  at  discrete  time 
steps,  using  the  full  MCELL  package,  is  useful  and  results 
are  illustrated  in  section  1'. 

The  TERRYS  radiation- ioni cat i on  package  utilizes 
R.  Terry'-  •  curve  fits  to  D.  Boston's  MCELL  computed 
optically  thin  volume  emission  rates  for  aluminum  or 
argon . ^ 0 ’ b ^  The  radiation  is  divided  into  the  categories 
of  line  L- shell  with  he  <  1  keY,  line  K- shell  with  hv  N  1  keY, 
free  bound  continuum  and  free-free  Bremsstrahlung  continuum 
radiation.  The  optically  thin  line  and  free-bound  emissions 
are  fitted  to  the  following  functional  relation: 


l. 
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are  shcrt  and  icr.g  path  optical  aepths 


described  below.  The  cptical.lv  thin  -  creouency  emission 
coefficient  i s 
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and  1. •  are  entered  as  a  table  for  1C 


1° 

10x-  H: 


The  i-  frecuenc>'  optical  depth,  :or  either  p2th 

■  ^  h 

length,  is  given  by,  for  emission  at  the  i *"  cell, 
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is  the  absorption  coefficient.  In  performing  the  t ^ 
integration  the  variability  of  7g ,  ne,  nj,  and  3j  is 
accounted  fer  along  the  path  lengths  by  employment  of  the 
previous  time  step's  solution  parameters. 

The  short  and  long  optical  depth  path  integrations 
are  carried  out  from  the  individual  cell  in  two  opposing 
directions  to  the  boundary  of  the  cylindrical  plasma.  The 
same  idea  in  estimating  opacity  is  also  used  for  line  and 
bound- free  continuum  transport: 


4  7 


for  either  path  length.  The  I'cppler  broadened  cross 
section  is  used  for  line  transport: 


o  _  a 


For  K-sheli  radiation  it  is  assumed  that  the  wavelength 
is  evaluated  at  a  single  value  of  1.0  keV  photons 
and  for  L-sheii  radiation  it  :s  assumed  that  0.5  keV  single 
wavelength  emission  is  appropriate.  Bound-free  opacity 
is  usually  small  for  the  conditions  of  interest;  the  cross 
section  used  is 
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where  £(!„';  is  a  pol vnomial  function  of  the  ion  number: 

X 

it  has  been  determined  by  a  fit  to  multi-cell  input  data. 

The  opacity  of  a  cylindrical  plasma  affects  the 
resulting  net  radiant  emission  by  removal  of  photons  prior 
to  their  escape  through  the  plasma  boundary.  Probability  of 
escape  functions  of  optical  depth  are  f i ts  to  soluti ons  of 
the  radiation  transport  equation  ’  ^  and  are  of  good 
accuracv  for  cvlindrical  transport  over  a  wide  range  of 
The  total  loss  from  all  radiant  moJes  at  each  cell  is  then 
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where  P(t'-  are  the  probability  of  escape  functions  defined  for  bound- 
free  radiation  bv  the  f  c-1 1  ovine  '  iron.  R .  Terrv,  private  communi  ca  t :  on 


.  ?(:;■ 


for  ,S?  -  5.0  . 


0  2 


p  r  *■ 

P f  .  =  '1.0  -  8  .  2  o  8  ,\  10  ‘  t  ,  ;  '  i  .  0  f  0.426  t 
S  5,5. 

for  ~  ^  <  5,0  . 

For  line  radiation  the  Doppler  profile  results  of  Ref.  S 
(Appendix  A;  are  utilised  by  R,  Terry  in  the  form 


=  1 . 0 1  |i . 0  - 
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The  ionization  model  utilizes  the  near-coronal 
model  Shearer- Earnes 1 ^  expression  (Tg  in  eV) 


“I 


2b.  0 


T; 1000.0 
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.14  ! 


for  the  ionization  number,  where  1.4  is  the  element  atomic 
number,  while  the  internal  electron  energy,  given  by 


E  =  4-  n  k  T  -  ET-n  . 
e  e  e  Tup 


CD 


utilizes  a  hydrogenic  assumption  for  the  total  chemical 
potential  Ey^p,  i.e.. 


"I  'I 

cTCP  nI  5Z  !F(‘i  '  ~  T"‘  J 
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where  R.  =  13.0  el'.  The  ionization  level,  given  by  Eq .  (34  , 
is  accurate  up  to  n^  =  10"^  cm.  but  should  be  close  enough 
for  densities  perhaps  a  factor  of  3  higher  than  this, 
covering  the  range  of  collapse  conditions  found  in  Black¬ 
jack  III  or  V  shots.  The  y  adiabatic  exponent  is  defined 
in  terms  of  the  electron  thermal  plus  ionization  energy  by 
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As  stated  in  section  II  the  adiabatic  exponent  is  evaluated 
from  the  temperature  and  ionicaticn  number  locally  at  each 
spatial  cell  and  timestep  for  utilization  in  the  conserva¬ 
tion  of  electron  temperature  and  total  internal  energy 
equations;  the  possible  influence  of  spatial  and  temporal 
gradients  of  y  evaluated  from  the  respective  T&  and 
gradients  is  not  taken  into  account. 

Most  of  the  computations  illustrated  in  this 
report  assume  microscopic  classical  transport  processes 
although  some  consideration  is  given  to  anomalous  transport, 
especially  resistivity.  The  classical  resistivity  is^10-' 

ta  =  1.15  x  10'14  ‘  T "  '  ~  i-L  (seel  (381 

1  e  •  - 

(  -)  •; 

uhe  re  '  * 

i:J.  =  23.  .1  -  1.15  log2p(ne!  -  3.45  1  Ogj  p  (T& ) 

for  T  s  50  eV  , 
e 

=  -5.3  -  1.15  log1()(ne;  *  2.3  logJp(Te) 

for  T  >  50  eV  . 
e 
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an-  i  i  ■  -i  ,  i  a  0  i 

w.iere  *  T  is  the  classical  value 
e  I 

=  3.44  x  102  Tg/:/(r.  Zj  i7J\ )  (sec)  (40a 

if  the  ion  sound  speed 


C, 


c  a  i  •  hh 
V  mic' 


(40b 


is  greater  than  the  absolute  value  of  the  drift  speed  U-, 
given  by  Eq.  (10),  and  <  T^.  However,  if  the  opposite 
conditions  prevail  at  a  local  radial  cell  the  rej  is 
calculated  from 


l  JOc 


where  u,p  ,  the  ion  plasma  frequency,  is 
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Although  this  formulation  is  superior  to  the  previous 
resisitivty  multiplier  scheme  the  inability  of  the  code 


tc  yield  verv  realistic  ion  temperatures  independent  cf 
the  electron  values  has  prevented  its  use  for  more  than  a 
few  test  cases  (Section  IV). 
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Section  IV 


BENCHMARK  CALCULATIONS 

Included  in.  this  section  are  ideal  gas  cylindrical 
diaphragm  implosion  as  well  as  magnetic  field  diffusion 
benchmark  results.  In  addition,  some  comparative  implo¬ 
sion  results,  with  and  without  radiative  coupling  or 
reconing,  are  illustrated.  It  was  decided  to  perform 
these  calculations  using  interesting  subsets  of  the  complete 
model  so  as  to  strengthen  the  level  of  confidence  in  the 
full  code.  The  first  illustrated  benchmark  is  the  hydro¬ 
dynamic-only  subset  where  the  voltage  discharge  and  induced 
magnetic  field  contributions  are  removed.  The  driver  of 
the  fluid  dynamics  is  t=0  radial  distributions  of  density 
and  pressure.  Implosion  is  generated  by  the  so-called 
cylindrical  diaphragm  problem  with  an  assumed  step  density 
and  pressure  rise  of  a  factor  of  four  outside  of  a  half-maximum 
radial  point. 

Criteria  to  observe  within  the  implosion  are  whether 
temperature,  pressure,  and  density  iumps  across  the  inwardly 
propagating  shock  and  contact  d i scont i nu i t i es  are  correctly 
transported  by  the  FCT  algorithm;  these  iumps  can  be 
calculated  from  steady  planar  one  -  dimens iona 1  theory 
coordinate  propagating  with  respect  to  the  discontinuity! 
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before  the  advancing  shock  gets  reasonably  close  to  the 
cylindrical  centerline.  Beyond  that  comparison  has  to  be 
made  Kith  established  numerical  solutions  in  the  literature, 
such  as'Lapidus.  The  solved  equation  subset  comprises 

Eq.  (1),  Eq.  (2)  without  the  iB/c  term,  Eq.  (5)  without 
the  i e  and  terms,  the  single  specie  internal  energy 

balance 


and  p  =  cRT,  for  the  five  unknown  c,  v,  E,  p  and  T. 

Illustrated  in  Figures  2  through  3  are  normalised  radial 
pressure  and  density  distributions  evaluated  at  four  dif¬ 
ferent  normalised  times.  Normalisation,  for  density  and 
pressure,  (s*,p*)  is  to  the  respective  t=0,  r=0  magnitudes, 
while  for  time  (t*)  it  is  to  the  unit  diaphragm  radius 
divided  by  the  sound  speed  at  t=r=0.  Included  in  Figures  3 
and  5  are  the  artifical  viscosity,  y  =  ].4  Lax-h'endroff 
finite  difference  solutions  obtained,  with  a  Cartesian  x-y 
mesh,  by  Lapidus.  Correlation  is  reasonably  good  within 

the  inner  range  of  the  implosion;  as  the  FCT  code  employs 
an  outer  boundary  condition  of  3/crisv}  =  0  while  the 
Ref.  12  solution  considers  pv  =  0  here,  considerable  discrep¬ 
ancy  in  the  outer  region  is  understandable.  The  ;ero  gradient 
boundary  condition  is  more  correct  for  r-pinch  implosions. 
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V.'hat  is  of  most  interest  is  the  shape  and 
inward  progression  of  the  pressure  2nd  density  shock  and 
contact  discontinuity  fronts;  the  FCT  calculated  results 
are  slightly  more  spread  out  than  the  Ref.  12  solutions. 

The  difference  is  largely  due  to  the  factor  of  five  decrease 
in  ir  increment  sice;  the  relatively  coarse  mesh  is  still 
adequate  to  illustrate  the  disturbance  propagation.  The 
FCT  implosion  is  also  slightly  slower  in  running  into  the 
centerline  than  the  Lapiaus  model  results:  the  outer  cere 
gradient  boundary  condition  on  ov,  leading  to  some  mass 
blowoff,  would  retard  the  inward  progression.  It  seems 
as  though  the  Ref.  12  pressure  and  density  results  do  not 
line  up  as  well  as  they  could  and  that  the  Lapidus  density 
calculations  are  running  a  bit  ahead  of  the  corresponding 
pressure  profiles.  The  FCT  density  profiles  may  therefore 
be  more  accurate  even  with  the  reduced  number  of  computation 
points . 

For  planar  flow  initiated  by  the  bursting  of  a 
diaphragm,  where  subscript  1  denotes  the  region  ahead  of 
the  propagation  shock,  subscript  2  the  region  immediately 
behind  the  shock,  subscript  3  the  region  behind  a  contact 
discontinuity  and  subscript  4  the  ambient  region  behind  an 
isentropic  expansion  fan,  the  Pj/P}  ratio  is  given  by^-5-1 


which  is  obtained  from  the  facts  that  pressure  is  uniform 
across  the  contact  surface  (p,  =  p - ) ,  there  is  an  isentropic 
change  between  regions  5  and  4,  and  from  the  initial 

conditions.  Finally,  for  the  perfect  gas 
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t  p5/cI  =  (P;/P1) (T1/T.)  .  (42a) 

Included  in  the  bottom  row  of  Table  1  are  the  numerical 
,  magnitudes  obtained  from  Eqs.  (42)  through  (42a).  The  other 

rows  of  Table  1  collect  results  at  four  t*  normaliced  times, 
given  by  the  numerical  solution  to  the  conservation  equation 
Any  discrepancies  with  the  above  results  are  due  to 
cylindrical  coordinate  geometric  focusing  and 
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the  likelihood  that  a  mesh  with  only  39  points  is  slightly 
coarse.  The  numerical  results  are,  as  would  be  expected, 
higher  than  the  planar  solution  but  are  close  enough  to 
indicate  a  suitable  numerical  hydro-only  algorithm. 

The  next  subset  of  the  program  to  benchmark  is  the 
diffusion  of  acimuthal  B-field,  described  by  Eq .  (13). 
Hydrodynamic  plus  B-field  solutions  will  first  explore  the 
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diffusior.less  limit  cf  r  =  £  =  0,  i.e.,  bcth  infinite 
conductivity  and  u>  -reT.  The  resulting  magnetic  flux 
equation , 


c  5 
ct 


(43 


is  of  similar  form  to  the  planar  mass  conservation  expression 
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Numerical  solution  comparisons  will  only  be  made  with 
analytical  estimates  across  a  plane -propag2t ing  wavefront. 
The  shock  relation  corresponding  to  Eq.  (45)  is 


Blvl  =  B;v: 


(45) 


while  mass  conservation  through  the  shock  surface  is  given  by 


clvl 


(46) 


Combination  of  the  above  two  relations  gives 


fc-i  =  B^fco/Cj) 


(471 


Likewise,  letting  E  =  ac,  where  a  is  to  be  determined,  we 
find  that  Eq.  (45)  becomes,  with  employment  of  Eq.  (44), 
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Therefore,  a  =  5/c  should  be  constant  in  the  di f fus i onles s 
limit  Kith  respect  to  a  coordinate  system  fixed  at  the 
shock.  The  wave  propagation  solutions  will  be  the  same  as 
the  previous  cylindrical  diaphragm  model  except  that,  at 
t  =  0,  a  radially  uniform  10''’  Gauss  magnetic  field  is  imposed 
(both  within  and  outside  the  diaphragm  location).  It  is 
first  of  interest  to  see  whether  c  and  the  E  field  tangen¬ 
tial  to  the  wave  surface  are  numerically  propagated  in  a 
similar  manner,  in  particular  behind  the  initial  shock 
front.  Due  to  the  rapidity  of  B  field  propagation  induced 
by  the  electrical  resistivity  and  temperature  gradient 
contributions  Eq.  (15)  is  solved  implicitly  while  the 
cylindrical  version  of  Eq .  (44)  is  solved  by  explicit  FCT. 
The  di f fusionless  limit,  Eq.  (43',  is  also  solved  implicitly 
so  some  discrepancy  may  enter  here. 

Plotted  in  Fig.  b  is  a  =  B/c  which,  from  Eq.  (48), 
should  be  constant  in  the  diffusioniess  limit  with  respect 
to  a  coordinate  system  fixed  at  the  shock.  As  is  seen  the 
variability  of  a  over  much  of  the  time  scale  is  not  too 
large  so  that  implicit  B-propagation  is  adequate.  The 
similarity  breaks  down  at  longer  times,  as  the  front  nears 
the  axis;  numerically  the  cylindrical  form  of  mass  conserva¬ 
tion  is  employed  rather  than  Eq.  (44  .  From  Eq.  (4~) 

B-n/c-,  =  3.15  x  10'  Gauss  cm~/gm  for  all  times  so  that  the 
cylindrical  numerical  solution  averages  about  Si  beneath 
this  value  up  to  1.5  x  10  ^ 


♦ 


ol 


sec . 


A  sign!  f leant lv  different  temporal  post-shock 
-field  distribution  occurs  when  classical  magnitudes  c 
esistivity  and  thermal  flux  influence  are  incorporated 
..s  shc;wn  in  Table  1,  the  E  magnitude  rapidly  drops,  wit 
some  degree  of  oscillation,  due  to  the  small  time  scale 


he  diffusion  process. 


the  near  equal  time  pressure  profiles  are  relatively  small 
except  ahead  of  the  pulse  leading  edge;  here  there  is  a 
precursor  pressure  perturbation  that  grows  slightly  with  time 
The  next  set  of  benchmark  computations  involves  the 
combination  of  the  full  B-field  equation,  Eq.  ( 1 5 j  ,  and  the 
mass,  momentum,  internal  energy  and  electron  temperature 
equations,  Eqs.  si.  .  U.'  2n-  l*4  •  ^  constant 

„  _  2  f, *  adiabatic  exponent  is  assumec  with  nc  radiation 


transport;  nowever, 


there  is  an  induced  1  : ;e;c  resulting 


from  a  discharge-generated  hydrodynamic  collapse.  n&  anc 

are  calculated  from.  Eqs,  lilt  anc  i. -■>•*  .  respec  >.ii  el.»  . 

"  I 

The  idealized  voltage  —  time  shape  illustrated  in  rig.  8  is 
employed,  similar  to  the  one  used  in  the  Ref.  11  .simula¬ 
tions,  and  the  the  three  resistivity  models  — classical  , 
corona  multiplier  {Ecs.  (.59;  and  (.59a)  with  A  =  49  and 
E  =  10*“},  and  the  Sagdeev  instability  form  {Eqs.  (40), 

(40c) }  are  utilized. 

The  initial  conditions  are  radially  uniform  electron 
and  ion  temperatures  of  1  eV,  and  an  initial  annular  Gaussian 
distribution  of  ion  number  density 


r‘I  =  nI„ 


exp  \  -  ( r - r  . " / R ' 


1  C  -  C*  j 

where  the  assumed  conditions  are  nT  =10  cm  and 

1  >  II]  d  A 

5  x  1018  cm"'\  r  =  0.85  cm  and  R  =  0.25  cm.  The  discharge 
c 


length  :5  5.  i  cm  ar.J  the-  :  l^sta  ;5  assumed  to  be  A r  with  •» 
remaining  at  l.tr".  lie  resistance  and  inductance  of  the 

*  *S 

external  circuit  are  assum.eu  tc  re  '  .  S  chr:  and  1.5  x  1C  he 
respectively,  .a  lues  characte. . « c ic  cl  the  Pi  then  machine 
at  Physics  International. 

Figure  9  illustrates  the  maximum  racial ly- dependent 
plasma  density  cut  tc  a  maximum  tine  c:  ICC  nsec  for  the 
three  resistivitv  models  and  two  maximum  r. .  values.  Kith 

i 

1  V  .  • 

r.  -  „  =  r  x  lr-  "■  cm  '  tr.e  initial  mass  ner  unit  length  is 

*  ’ max 

p.PS  x  10  "  cm  cn  while  with  n-  =  10*g  cm.  *'*  the  mass 

*  -  max 

per  unit  length  is  1.59  x  IP  "  cm.  cm.,  which  is  comparable 

tc  the  Ref.  11  value  and  of  an  appropriate  magnitude  for 

a  maximum  "tailored"  implosion  for  these  machine ■ paramet ers . 

The  smaller  mass  yields  an  im.plcsicn  within  the  time  span  of 

the  voltage  pulse  while  a  tendency  towards  collapse  for  the 

larger  mass  only  occurs  after  pulse  cessation.  The  ion 

temperature,  calculated  from  Eq.  (5.:,  is,  for  the  most  part 

near  or  somewhat  larger  than  T  ,  so  that  the  Sagdeev  model 

e 

yields  values  usually  of  the  classical  magnitude;  however, 
the  density  peak  is  close  to  the  50  times  classical  corona 
multiplier  model.  This  anomalous  resistivity  multiplier 
model  yields  an  implosion  as  dynamic  as  the  classical  r 
solution;  indeed,  for  the  smaller  mass  case  the  r,  implo¬ 
sion  is  1"  nsec  more  rapid  although  the  peak  density  is 
somewhat  reduced.  This  result  is,  however,  in  contrast  to 
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the  softer  implosions  produced  with  resistivity  multiplier 
solutions  when  radiation  transport  is  included  (Ref.  11 
and  Section  v). 

Illustrated  in  Fig.  10  are  the  radial  distributions 
of  Tg  for  the  lower  mass  distribution  at  the  times  of 
maximum  compression  (depending  on.y).  The  temperature 
peaks  in  this  radiationless  calculation  are  immediately 
outside  the  region  of  rapid  density  fall-off,  as  seen  by 
comparison  of  Figs.  10  and  11.  Maximum  attained  densities 
are  not  as  high  or  the  width  of  the  collapsed  tone  not  as 
narrow  as  would  be  found  m  a  coupled  radiation  emission 
solution,  however.  In  contrast,  the  higher  mass  level 
(Fig.  12)  weakly  implodes  the  discharge  resulting-  in  a 
relatively  ineffective  collapse.  It  is  of  interest  to  note 
that  the  lack  of  significant  collapse  enables  the  202  nsec 
time  span  to  be  covered  with  500  minimum  Courant  condition 
time  steps  while  a  maximum  of  160  nsec  is  spanned  with 
950  Courant  condition  steps  when  a  collapse  is  achie\ed. 

It  is  noted  that  the  inclusion  of  radiation  and  the 
generally  more  restrictive  internal  energy  change  timestep, 
Eq.  (25d.l,  results  in  roughly  a  tripling  of  the  timestep 
number  (Section  V). 

Referring  again  to  Fig.  9  we  see  that  the  implosion 
time  for  the  five  times  heavier  mass  is  about  40$-  greater 
than  the  implosion  time  for  the  lighter  mass  (measured  from 
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the  beginning  of  collapse).  It  is  of  interest  to  see 
whether  a  simple  dynamic  pinch  model,  suggested  by  A.  Monoeli 
yields  anything  like  this  value.  If  we  assume  that  the 
discharge  acts  through  a  cylindrical  plasma  annulus  of 
thickness  -  ,  length  L,  and  radius  r,  such  that  p  =  (M/L)/(2r:r 
J  =  I  /  2  -  r :  and  B  =  2 1 / c r ,  then 


d-r  J  x  E 


(50) 


yields,  upon  collapse  from  r  =  r  to  r  =  0  over  a  timespan  t 


t  =  r  JJ\  cV^ 

\  L  21  • 


(  5  1  ; 


The  implosion  time  ratio  is  then  proportional  to  the 
square  root  of  the  mass  ratio,  within  this  simple  model. 

As  the  actual  correlation  is  closer  to  the  0.25  power  of 
the  mass  ratio,  simple  modeling  of  these  implosions  is 
of  more  qualitative  rather  than  quantitative  value. 

Test  calculations  evaluating  the  suitability  of 
dynamic  retoning  are  performed  using  the  Physics  Interna¬ 
tional  Pithon  Shot  ?1400  Argon  gas  puff.  The  TERRYS  two- 
path  radiation  package  is  employed,  contrary  to  the  employ¬ 
ment  of  the  ATPHVS  routine  on  the  same  data  last  year. 

The  argon  fits  to  the  optically  thin  radiation  parameters 
are  as  follows,  from  R.  Terry  (private  communication'. 


bb 


In  Eq .  (.50! 


X 


0.60,  M  =  14  for  L- she  11  line  radiation 
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The  t  =  0  ,  5  cm  long  initial  gas  puff  ion  number  densitv, 
over  the  59  cells  of  the  computation,  is  illustrated  in 
Fig.  13.  The  dashed  line  near  the  outer  edge  is  an  extra¬ 
polation  past  the  last  measured  data  point.  Input  is  bv 
simple  specification  of  c  *»  nx  nij  at  the  39  cell  midpoints. 
The  radially  uniform  temperature  at  t=0,  for  both  electrons 
and  ions,  is  1.0  eV. 

The  first  presented  calculations  assume  the  spatiallv 
uniform  cell  case.  Plotted,  at  the  two  times  of  81."  and 
101.:  nsec,  are  nj,  ng,  Tg  and  normalised  Fp,^  (total)  as 
a  function  of  radius.  Also  plotted,  as  a  function  of  time, 
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is  the  total  normalized  radiated  energy  per  unit  discharge 
length  out  to  t  =  125."  nsec.  Figure  14  illustrates  the 
number  density  distributions  at  both  times.  The  peak  of 
the  number  density  distributions  has  not  moved  appreciably 
inward  at  81.67  nsec  but  the  magnitudes  at  and  within  the 
peak  radii  have  increased  somewhat  at  the  expense  of  the 
outer  zones.  The  first  complete  implosion,  defined  at  the 
point  beyond  which  the  net  velocity  begins  to  reverse  itself 
near  the  discharge  centerline,  occurs  at  near  101.2  nsec. 

It  is  seen  that  steep  centerline  density  peaks,  extending 
only  a  few  cells  in  width,  have  formed  at  this  point. 
Actually  these  sudden  centerline  peaks  have  formed  only 
over  the  previous  3-4  nsec  of  implosion  time.  The  density 
steps  are  largely  due  to  the  retention  of  perturbations 
which  is  a  feature  of  the  FCT  algorithm,  as  well  as  a 
slightly  insufficient  cell  number.  The  relatively  flat 
floors  of  density  are  essentially  imposed  at  a  magnitude 
that  is  5  x  10  0  of  the  average  value  over  the  total  dis¬ 
charge  width.  A  comparatively  high  centerline  density 
(within  an  order  of  magnitude  of  the  illustrated  values) 
is  maintained  for  at  least  25  nsec  beyond  the  101.2  nsec 
point;  however,  the  radial  velocity  distribution  becomes 
quite  discontinuous  with  probable  formation  of  hydrodynamic 
instabi 1 ities . 
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Figure  15  shows  the  electron  temperature  distribu¬ 
tions  at  the  same  times.  It  is  seen  that  the  outer  tones 
of  the  implosion  heat  up  to  more  than  10  keV  coronal  tempera¬ 
tures;'  densities  are,  for  the  most  part,  correspondingly 
low  at  these  positions.  Near  the  maximum  compression  the 
high  density  centerline  material  heats  up  to  over  1  keV’ 
temperatures,  yielding  the  maximum  radiative  output. 

Figure  16  illustrates  the  radial  normalised  total 
radiative  power  distributions  at  the  same  times.  The 
radiative  peak,  at  81.6/  nsec,  is  at  near  the  maximum  re¬ 
location,  with  a  T  value  of  about  400  eV.  The  maximum 

e 

radiative  location  is  the  centerline  at  the  101.2  nsec 


initial  compression  time.  Figure  17  shows  the  to.tal 
normalised  radiant  energy  per  unit  length  for  both  the 
uniform  cell  width  case  and  a  resoning  results,  discussed 
later.  It  is  seen  that  increases  at  a  relatively  slow 

rate  after  the  first  collapse,  although  running  out  several 
hundred  more  steps  (it  took  629  steps  to  reach  81.6“  nsec, 
1119  steps  to  get  to  101.2  nsec,  and  2000  steps  to  reach 
125.7  nsec}  beyond  the  125. 7  nsec  stop  may  attain  a  measur¬ 
able  increase.  However,  as  indicated  before,  dis¬ 

continuities  in  the  radial  velocity  profile  point  to  an 
unstable  compressed  region,  with  a  1-D  modeling  being 
questionable. 
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Illustrated  in  Figs.  17  through  21  are  the  results 
of  a  retoned  calculation  with  M°  =  0.025,  b  =  0.03,  c  =  0.05 
{in  Eq.  (29)}.  The  solution  remains  stable  with  a  maximum 
divergence  from  energy  balance  of  2 %  (after  900  time  steps 
the  error  is  significantly  smaller).  Going  first  to 
Fig.  IS  we  see  the  cell  width  distributions  at  t  *  82.01 
and  100. "  nsec.  These  times  are  comparable  to  those  chosen 
in  the  uniform  increment  solution,  but  it  is  interesting  to 
note  the  different  number  of  time  steps  required  to  reach 
these  times  (1124  and  1649  versus  629  and  1119  for  the 
uniform  cells).  Comparison  of  the  lt  distributions  in 
Fig.  IS  with  the  associated  nj  distributions  in  Fig.  19 
show  that  appropriate  tracking  (smallest  cells  with  largest 
density  and  vice  versa)  is  produced.  The  average  cell 
width  is  calculated  from  Eq .  (27a).  At  the  same  implosion 
times  the  computed  outer  radius  of  the  discharge  is  about 
10%  smaller  in  the  retoned  solution  compared  with  the 
uniform  cell  model. 

Comparison  of  the  Figs.  14  and  Fig.  19  number 
density  shape  distributions  show  general  similarity  but 
sufficient  differences  to  indicate  some  cell  width  effects. 
In  particular,  the  retoned  maximum  centerline  n-j  at  the 
collapse  point  is  about  SOt,  higher  than  the  uniform  cell 
result.  Curiously  enough,  as  the  retoned  centerline 
temperature  is  somewhat  lower  the  radiated  power  at  this 
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point  is  about  the  same,  irrespective  of  the  model  (Figs, 
and  21).  The  narrowness  of  the  maximum  density  distri¬ 
bution  near  the  center  is  not  really  reduced  with  reconing 
although  it  should  be  noted  that  the  extreme  density 
"spike"  is  largely  formed  in  the  last  2  nsec  before  the 
collpase;  prior  to  this  the  density  and  radiation 
distributions  are  somewhat  less  peaked  near  the  center. 

Returning  to  Fig.  i"  we  see  that  the  radiated 
energy  output  obtained  with  the  reconed  mesh  is  comparable 
with  the  previous  result;  however,  at  times  beyond  the 
collapse  the  reconed  magnitudes  can  be  up  to  20fc  below  tha 
obtained  with  the  uniform  mesh.  In  summary f  comparison  of 
Figs.  14  and  19  for  number  density,  Figs.  15  and-  20  for 
electron  temperature ,  16  and  21  for  radiative  power  and 
the  separate  energy  distributions  on  Fig.  17  illustrates 
the  influence  of  a  moderate  reconing  algorithm. 

In  conclusion,  the  varous  subsets  of  the 
code-cylindrical  implosion  dynamics,  magnetic  flux 
generation,  dynamic  cell  reconing,  nonequilibrium  radia¬ 
tion  coupling,  and  resistivity  transport  modeling  have 
been  checked  out  separately  and  in  some  combination;  it 
is  felt  that  full  modeling  with  realistic  input  data  can 
give  suitable  theoretical -experimental  correlation. 


Section  V 

SAMPLE  CALCULATIONS 

The  coupled  hydrodynamic - radiat ion  calculation? 
are  illustrated  for  two  separate  situations.  The  first 
is  a  comparison  with  the  SPLAT  code  implosion  results 
due  tc  R.  Terry1"51,  following  the  original  analysis  of 
Tidman  " ‘5"} :  in  this  code  the  plasma  produced  in  the  dis¬ 
charge  through  an  A I  wire  cage  is  modeled  by  a  radially 
isothermal  analysis  with  radially  integrated  representa¬ 
tions  of  mass,  momentum,  and  energy  across  an  assumed  core 
and  corona  toning.  This  is  a  so-called  zero-dimension 
transient  analysis  in  contrast  to  the  one -dimensional 
solutions  illustrated  here.  The  second  comparison  set  is 
a  numerical  representation  of  a  Maxwell  Laboratory  Black¬ 
jack  Y  experimental  con  figuration;  a  variation  of  plasma  mass, 
for  both  At  and  Ar,  is  employed  so  as  to  compare  the  space 
and  time  integrated  radiation  yields.  In  addition,  radial 
distributions  of  temperature,  radiative  power,  and  density 
are  illustrated  for  near  the  optimum  mass  at  the  pinch 
collapse  time. 

Figure  22  shows  the  voltage-time  profile  employed 
in  the  SPLAT  calculations.  The  maximum  voltage  is 
4. 5 x 10^ Y,  attained  in  150  nsec.  The  quantities  would  be 
experimentally  maximum,  and  indicate  an  optimistically 
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energetic  profile.  External  circuit  parameters  are  a 

.  o 

0 .  5  7  generator  impedance  Ic  and  a  13.5  x  10  "  h  diode 

inductance.  Although  a  variety  of  original  wire  cage 

radii  were  investigated  in  Ref.  4,  ranging  from  R"  =  0.7 

tc  1.6  cm,  only  Rc  =  1.0  cm  will  be  explored  here.  It  is 

assumed  in  all  the  wire  calculations  that  there  is  an 

acimuthai  symmetry  in  the  plasma  prior  to  an 

appreciable  implosive  collapse.  The  calculations  assume 

three  different  masses,  made  up  of  6,  12,  and  24  wires  of 

1  mil  diameter  and  3  cm  length.  The  mass/wire  of  Ak  wire 

is  4.10  x  10  D  gm  sc  that  the  three  separate  masses 

analyzed  are  2.46  x  10  4,  4.92  x  10  4,  and  9.S5  x  10  ^  gm , 

respectively.  It  is  assumed  that  a  39  cell  numerical 

analysis  is  adequate  and  that  most  of  the  initial  plasma 

density  is  included  into  the  two  cells  nearest  r  =  1  cm, 

the  total  computation  radial  extent  being  near  2  cm  at 

t  =  0.  The  initial  six  wire  plasma  density  distribution 

is  then  1.0  x  10*  gm/crn-5  for  .0256  s  rc  ^  0 .  S 9  “ 7  cm, 

1.28  x  10  4  gm/cnT  for  .949  <;  r^  s  1.000  cm,  and 

1.0  x  10  '  gm/cm0  for  1.0516  s  r_  s  1 . 9 “ 5  cm,  r  denoting 

^  c 

the  cell  center  location.  Twelve  and  24  wire  cages  are 
corresponding  multiples  of  the  .949  s  r_  i  1.000  cm  density 
value.  A  small  amount  of  mass  is  therefore  assumed  to  have 
diffused  over  the  entire  0  <;  r  <;  2.0  cm  computational 
extent  at  the  initialization  of  calculations. 
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It  is  assumed  that  uniform  cell  width,  nc  density- 
dependent  reconing,  is  adequate  for  these  computations. 

The  initial  radially  uniform  temperature  is  10  eV. 
Illustrated  in  Fig.  25  is  the  temporal  distribution  of 
ncvmalited  energy  for  the  three  wire  numbers;  computations 
were  stopped  when  excessive  velocity  discontinuities  occur 
and  energy  balance  begins  to  break  down.  The  same  energy 
magnitude  normaliration  is  used  for  this  case  as  is  employed 
in  the  subsequent  Blackjack  V  representations. 

Table  4  shows  the  time  to  collapse  (pulse  time) 
for  the  three  masses  as  calculated  by  the  R.  Terry  SPLAT 
code  and  the  WHY RAD  results;  the  KHYRAP  computed  time  is 
the  location  of  the  first  peak  of  cross- sectionally  inte¬ 
grated  radiative  power,  simultaneous  to  the  maximum  yield 
located  near  the  centerline.  Also  included  in  Table  4 
are  the  ratios  of  total  radiated  energy  at  these  times 
divided  by  the  six  wire  value  for  the  KHYRAP  calculations, 
as  well  as  the  ratios  of  the  K-shell  line  radiation 
magnitudes  for  both  SPLAT  and  KHYRAP  computations. 


Table  4 

SPLAT  KHYRAP  K- Shell  K- She 11 


Wire  No. 

Pulse 

(nsec) 

Time 

it 

Pulse  Time 
(nsec)  ct 

KHYRAP 

E  ’E^ 

CRAP  CRAP 

KHYRAP 

erap  erap 

SPLAT 

F  /RP 

CRAP'  CRAD 

6 

98.9 

111.6 

1.0 

1.0 

1.0 

15.8 

10.2 

1  - 

114 .  “ 

8.4 

i:i'8  19.0 

1.6S8 

1.415 

2.02 

24 

123.1 

140.8 

4.196 

1.910 

2.6" 
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Direct  comparison  between  the  WHY RAD  anc  SPLAT 
results  could  be  better;  in  particular,  the  incremental 
time  change  between  successively  massive  collapses  is  nearer 
tc  a  linear  function  of  wire  number  with  WHY RAD  while  its 
a  sharply  decreasing  function  of  wire  number  with  SPLAT. 

The  h-shell  line  WHYRAD  yield  ratios  are  somewhat  lower 
than  the  comparable  SPLAT  results,  although  the  total 
radiation  yield  as  obtained  from  WHYRAD  steeply  increases 
for  11  wires  due  to  continuum  and  L-sheli  line  contributions. 
Nevertheless,  there  is  a  steady  vield  increase  with  both 
codes  as  the  wire  number  increases,  and  the  wire  cage  is 
not  yet  too  massive  sc  as  to  be  weakly  collapsible. 

Figure  21  shows  radial  electron  and  ion  number 
density  distributions  at  the  12  wire  collapse  point, 
t  =  125.8  nsec.  This  is  a  few  nsec  beyond  the  time 
selected  for  the  illustration  of  comparative  yield.  Fig¬ 
ure  25  illustrates  the  radial  distributions  of  normalited 
total  radiative  power  and  electron  temperature  at  the  same 
time.  The  temperature  near  the  centerline  is  relatively 
high  at  this  time;  the  local  yield  of  radiant  energy  subse¬ 
quently  cools  these  cells  to  under  100  eY  with  eventual 
reheating  due  to  thermal  conduction. 

The  Blackjack  Y  representations  illustrated  here 
have  been  presented  at  the  19S1  APS  Plasma  Physics  meet  ing .  ^ 3  ^ 
So  as  to  show  the  comparative  flexibility  of  the  WHYRAD 
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program  it  was  decided  tc  run  both  A;.  and  Ar  implosion 
examples,  using  the  same  voltage-time  profile  over  a  range 
ef  discharge  medium  masses.  The  Maxwell  Blackiack  V  machine 


cone itions ,  as 

applied  tc  an 

Ar  gas 

puff,  were  presented 

in  a  memo  by  Ric 

h  a  r  d  s  c  n  ,  e  t  a 

7  (16). 

*  *  J 

while  more  general 

machine  infermat 

ion  (external 

circui 

t  inductance,  impedanc 

is  contained  in 

an  R & D  Assoc  A 

ates  report  by  Fisher  and  See 

The  ini 

tial  radial  temperature  distribution  is 

assumed  to  be  a 

uniform  1  eV , 

wh.  lie 

the  timewise  diode 

voltage  distribu 

tion  is  given 

in  Table  5: 

Table  5 

t  (nsec 

V£  llO6  Vi 

t (nsec 

j  ve  UP6  V)  . 

0.0 

0.0 

85.0 

1 .  50 

42.0 

1.5 

:  95.0 

1.50 

50.0 

2.0 

112.0 

0.60 

60.0 

o  —  •- 

122.0 

0.40 

64.0 

2.80 

125.0 

0.00 

68.0 

,  -  - 
2  .  .  o 

140.0 

-0.  “5 

"5.0 

2.3“ 

150.0 

-1.00 

.  0 

2.00 

175.0 

-0.50 

80.0 

1 .  90 

200.0 

-0.10 

Linear  interpola 

tion  between 

these  points  accurately  model 

the  Blackiack  V 

voltage  pulse 

of  the 

Richardson ,  et  al  .  , 

experiments.  Figure  26  illustrates  both  the  voltage  and 
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current  distributions.  The  machine  impedance  and  indue- 

-  Q 

tance  are  0.58  f.  and  2  x  10  '  nenries,  respectively,  while 

-  0  -  Q 

the  diode  inductance  varies  from  15  x  10  '  to  2"  x  10  *  henries, 
depending  on  diode  geometry.  liayne  Clark  (private  communica¬ 
tion)  at  Maxwell  has  indicated  that  the  appropriate  diode 

-  o 

inductance  is  21  x  10  *  henries,  so  that  with  the  open 

circuit  prior  to  crowbarring  being  in  series, 

-  o 

L  =  25  x  10  "  henries. 

The  Ref.  16  initial  Ar  gas  puff  pressure  (cr  density) 
distributions  are  very  qualitative,  so  it  was  decided  to 
assume  a  general  Gaussian  relation  for  the  t  =  0  density 
distribution : 


c  =  :MAX  exp  {  -  E  (.  r  - 1 . 0  T  -  }  (52j 

where  1.0  cm  is  the  annular  radius  of  the  gas  puff  maximum 
density  in  the  Richardson,  et  al.  experiments.  The  density 
is  not  really  symmetric  about  this  point  as  more  material 
is  initally  present  throughout  the  outer  cones;  however, 
this  distribution  is  suitable  for  comparative  calculation 

v 

with  an  At  wire  cage  located  at  the  same  radius.  This, 
of  course,  also  assumes  that  the  wire  plasma,  upon  forma¬ 
tion,  collects  in  an  acimuthal  direction  prior  to  any 
appreciable  radial  collapse.  It  is  assumed  that,  for  the 

same  mass,  the  Ai  plasma  distribution  is  somewhat  narrower  \ 

around  the  Gaussian  peak  at  t  =  0  than  the  Ar  gas  puff  so 


,  S 
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that  =  5,  =  3.  Four  separate  ir.asses  are  employed 

for  At  and  a  subset  of  two  separate  masses  for  Ar ;  the 
four  masses  are 

M  =  50,  10,  150  and  300  ugm, 

the  lowest  mass  being  a  bit  small  for  wires  (=»  2  A£  wires)  , 
while  the  highest  mass  is  representative  of  a  12  wire  cage 
or  a  rather  dense  Ar  puff.  Kith  the  specification  of  the 
total  mass  then  in  Eq .  (51),  is  determined  by  numerical 

integration  of 

R° 

2r  [  e-BV-l.Or  rdr 
0 

where  R°  is  assumed  to  be  2.0  cm.  L  -  3.0  cm  is  the 
diode  spacing. 

The  flow  is  divided  into  45  ceils  out  to  a  total 
initial  radius  of  2  cm,  six  more  than  previously,  so  as 
to  hopefully  yield  a  bit  more  accuracy  near  the  maximum 
collapse.  In  this  regard  the  M°  power  in  the  density-cell 
width  retoning  algorithm  is  assumed  to  have  the  value 
0.025  so  as  to  obtain  a  50-40®  maximum  cell  width  differential 
(.probably  as  much  as  can  be  obtained  with  the  present  code'.. 
Classical  transport,  including  the  resistivity,  will  be 
assumed  in  the  initial  group  of  calculations. 
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Figure  2"  shows  aluminum  ion  and  electron  number 

density  radial  distributions  at  t  =  0  and  two  subsequent 

times  for  M/L  =  5  x  10  5  gm/cm  (M  =  150  ugm)  .  This  mass 

magnitude,  as  will  be  shown  later,  is  about  the  optimum 

implosion  value.  The  105.7  nsec  time  is  prior  to  the 

rapid  pinch  collapse,  with  resulting  energetic  radiation 

yield,  while  118.5  nsec  is  after  the  collapse.  The 

dynamic  sweeping  in  of  n^  between  these  two  times  is 

evident.  The  collapsed  high  density  region  diameter  of 

•>0  -  3 

about  1  mm  is  realistic  as  is  the  10"  cm  ion  density 
y  -  * 

and  10“  cm  "  electron  density.  Illustrated  in  Figure  28 
are  the  temperature  and  normalized  radiative  power  radial 
distributions  at  the  same  mass  level  and  implosion  , times  ; 
the  arrows  point  to  the  ordinate  associated  with  the 
particular  curve.  The  centerline  temperatures  heat  up  to 
near  1  keV  (from  the  initial  1  eV  uniform  temperature) 
while  low  density  corona  temperatures  are  >  10  keV.  The 
radiative  power  yield  after  collapse  is  strongly  peaked 
near  the  axis,  however,  as  the  density  is  so  low  in  the 
corona . 

Figures  29  and  30  show  the  time  dependencies  of 
the  cross-sectional  integrated  normalized  radiative  powers 
for  the  four  masses  of  interest  (M/L  =  1.66",  3.333,  5.0 
and  10  x  10  ^  gm/cm).  These  distributions  show  a  multiple 
peaking  behavior;  upon  compression  there  is  a  heating  up 


of  the  inner  cells  resulting  in  large  radiation  outputs 
with  subsequent  cooling.  The  radiative  output  in  this 
high  density  region  then  falls;  however,  thermal  diffusion 
processes  subsequently  induce  a  reheating  with  an  eventual 
repetition  of  the  radiative  peaking.  It  is  possible  that 
incorporation  of  viscous  momentum  transport  and  energy 
dissipation  would  have  seme  damping  effect  on  the  radia¬ 
tive  peaking.  It  is  seen  that  the  initial  collapse, 
illustrated  by  the  first  power  peak,  progresses  back  from 
90  nsec  to  140  nsec  as  the  mass  is  increased.  If  it  is 
assumed  that  an  optimum  tailored  implosion  is  at  the  maxi¬ 
mum  current  point  then,  from  Figure  26, the  initial  power 
peak  should  be  about  110-120  nsec.  The  M/L  =  5.0  x  10  5  gm/cm 
case  would  seem  to  be  the  most  optimum  although,  as  is 
illustrated  in  Figure  31,  the  normalised  radiant  energy  is 
slightly  maximised  with  the  3.53  x  10  5  gm/cm  mass  case. 

As  indicated  in  the  previous  computations  the  radial 
velocity  distribution  becomes  highly  di sc ont i nuous  in  the 
past-initial  collapse  phase  and  the  accuracy  of  the  hydro- 

dynamic  solution  rapidly  decreases.  They  dynamics  of  the 
1 

10  ■*  gm/cm  collapse  is  rather  weak;  probably  doubling 
this  value  would  prevent  any  collapse  at  all. 

The  maximum  computed  energy  yields,  as  given  in 
the  classified  DXA  final  report,  are  somewhat  smaller  than 
recent  experimental  results  (private  communication  from 
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J.  Peariman  and  K.  Clark);  however,  the  experimentally 
optimum  tailored  mass  seems  to  be  about  150  -_gm  or  a  little 
greater.  A  theoretical  yield  increase  could  be  obtained  by 
further  heating  of  the  center  two  or  three  cells  upon 
collapse;  relatively  modest  temperature  increases  would 
significantly  raise  the  radiation  level.  It  is  possible 
that  incorporation  of  viscous  dissipation  in  the  energy  equa¬ 
tions  would  also  help  in  this  regard,  and  it  is  also  possible 
that  evaluation  of  axially  dependent  hot  spots  by  two 
dimensional  modeling  would  yield  substantially  higher 
radiant  energy  magnitudes.  These  avenues  of  investigation 
should  be  pursued  in  future  work. 

Figures  52,  55  and  54  illustrate  area  integrated 
radiative  powers  and  radiative  energies  (both  normalised) 
for  A£  and  Ar ,  considering  the  lower  two  masses.  As  was 
mentioned  in  Section  IV  the  L-shelT  Ar  radiation  is,  at 
present,  incomplete;  therefore,  comparison  between  the  two 
elements  is  illustrated  by  removing  the  L-shell  contribu¬ 
tion  from  A£ .  It  is  interesting  to  note  that  the  temporal 
dependencies  of  the  integrated  radiative  powers  are  similar 
for  both  elements;  A£  has  a  higher  peak  at  collapse,  how¬ 
ever,  with  neglect  of  the  L-shell  radiation.  It  is  also 
interesting  to  note  the  collapse  times  for  the  two  AJL 
masses,  ~90  nsec  and  105  nsec,  are  the  same  as  the  pre¬ 
viously  illustrated  At  collapse  times,  including  the 
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L-s'nell  radiative  contribution.  Therefore,  any  reason¬ 
ably  modest  alteration  of  the  radiative  modeling  and  the 
resulting  yields  would  fortunately  not  have  too  much  of  an 
effect  on  the  relative  mass  tailoring  illustrated  here. 

The  multiple  peaking  behavior,  illustrated  in  the  previous 
power-time  distributions  (Figs.  29  and  30),  is  also  pre¬ 
sent,  in  particular,  for  A£  (Figs.  32  and  53).  The  actions 
of  thermal  diffusion  and  radiative  transport  are  an 
explanation.  Neglect  of  the  Ar  L-shell  radiation  seems 
to  produce  a  less  oscillatory  radiative  distribution  with 
a  more  rapid  post-initial  compression  drop. 

Figure  35  compares ,  for  an  M/L  =  5.353  x  10  ^  gm/cm 
A£  implosion,  the  effects  of  anomalous  resistivity  in  the 
normalized  radiant  energy.  Although  the  anomalous  resis¬ 
tivity  multiplier  model  initially  yields  a  higher  radiant 
energy  the  end  result  is  a  somewhat  softer  implosion  with 
a  501  reduction  in  total  yield. 

Figures  36  through  43  illustrate  L  and  K -shell 
aluminum  spectra,  line  plus  continuum,  as  obtained  from 
post-processing  at  the  two  separate  time  steps  shown  in 
Figures  27  and  28 ,  for  the  M/L  =  5.0  x  10* 5  gm/cm  mass  value. 
The  individual  temperature  and  density  distributions, 
over  a  sufficient  number  of  inner  cells,  were  inserted 
into  the  full  MCELL  multicell  coupled  radiation  and 
population  density  routine  by  D.  Duston.  The  R  value  is 
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the  radial  extent  of  the  computational  cells  inserted  into 
MCELL;  radial  extension  beyond  this  value  is  unnecessary 
due  to  the  optical  thinness  and  reduced  emission  of  the 
low  density  corona. 

The  L-shell  spectra,  prior  to  the  initial  collapse, 
is  shown  in  Figures  36  and  3"  over  the  photon  energy  range 
from  0.03  to  0.42  keV.  The  same  normalization  intensity 
value,  in  (U'/cnT/keV)  is  used  in  all  the  spectral  distri¬ 
butions.  Most  of  the  energetic  L-shell  lines  are 
Li-like;  it  is  seen  from  Figure  3"  that  a  normalized 
intensity  of  0.5  is  reached  with  the  2s-3p  Li-like  line  at 
0.256  keV  photon  energies.  The  2p-3d  line  is  only  a 
little  less  energetic.  Figure  38  illustrates  the  pre¬ 
collapse  less  complex  K-shell  spectra  made  up  of  AJc.*111 
H-like  and  AL*11  He- like  lines  as  well  as  the  continuum 
intensity  distributions.  The  most  intense  Lyman -a  line 
gets  up  to  a  normalized  intensity  value  of  0.6  at  1.728  keV 
photon  energy.  Figure  59  is  a  integrated  spectrum  that 
ideally  should  be  seen  by  X-ray  diodes;  for  example,  a 
diode  opaque  to  all  radiation  below  1  keV  would  present 
a  spectrum  having  the  illustrated  shape  for  E  >  1  keV. 

The  post  collapse  spectra  at  t  *  118.5  nsec  is 
illustrated  in  Figures  40  through  43;  a  number  of  individual 
line  intensities  are  markedly  increased  upon  collapse. 

As  an  example,  the  11-  like  2p-5d  line  now  emits  a  normalized 


intensity  of  about  15.0  at  hv  =  0.319S6  keV .  For  K-shell 

radiation  the  L  line  nuts  out  a  normalized  intensitv  of 
a 

about  120.0,  which  is  200  times  the  intensity  prior  to 

the  collapse.  A  similar  intensity  increase  of  about  two 

orders  of  magnitude  would  be  presented  by  XRD's,  as  is 

seen  when  comparing  Figures  39  and  43. 

Summarization  of  the  Blackjack  1'  -  simulated  numerica 

representations  should  include  the  following  statements: 

Centerline  ion  number  densities  implode  to  the 
*>n  -  3 

order  of  10“  cm  with  near  1  keV  temperatures  at  the 
same  location.  The  computed  high  density  region  is  around 
1  mm  in  diameter  upon  collapse. 

Collapse  times,  as  illustrated  by  temporal  varia¬ 
tions  of  the  radially  averaged  radiative  power  as  well  as 
the  radiative  energy  may  be  directly  compared  against  the 
peak  of  the  experimental  total  current  distribution. 

The  100  ygm  and  150  ygm  A 1  mass  implosions  present 
ootimum  tailoring  for  the  Blackjack  V  simulations.  The 
300  ugm  case  is  too  massive,  and  while  eventually  imploding 
at  =-  140  nsec,  the  yield  is  somewhat  reduced. 

The  energy  yields  exhibit  a  rapid  rise  of  3-5  nsec 
immediately  prior  to  collapse,  indicating  a  fast  creation 
of  the  centerline  density  peak.  Although  not  specifically 
illustrated,  the  computed  centerline  temperature  values 
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fluctuate  in  time,  falling  to  temporarily  low  values 
:  —  5 0  eV)  immediately  after  2  high  power  yield  peak,  with 
subsequent  reheating. 

Continuation  of  the  solution  to  times  much  beyond 
the  initial  collapse  is  of  reduced  practicality  due  to 
calculated  sharp  discontinuities  of  the  radial  velocity 
profile;  beam  instabilities  would  predominate. 

Comparison  of  continuum  plus  K-shell  Argon  and 
aluminum  implosions  indicate  fairly  close  temporal  power 
and  energy  distributions  with  Ai  eventually  attaining  a 
larger  magnitude  with  a  more  sudden  energy  jump  upon  collapse. 

Incorporation  of  a  multiplicative  anomalous  resis¬ 
tivity  in  the  coronal  region  of  the  plasma  yields  a  softer 
implosion  with  a  reduced  total  energy  output. 

Line  radiation  intensities  exhibit  a  rapid  increase 
upon  collapse  with  the  hydrogen- 1 ike  Lyman  lines  producing 
at  least  two  orders  of  magnitude  jumps  within  a  10  nsec  or 
so  time  span.  Post  processing  of  the  coupled  hydrodynamic - 
approximate  radiation  solutions  with  a  detailed  MCELL 
computation  furnishes  an  effective  comparison  diagnostic 
with  fast  time  scale  experimental  pulse  spectral 
measurements . 
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Section  VI 
CONCLUSION 

This  section  has  covered  the  basic  modeling  behind 
h'HYRAD  with  its  application  to  a  wide  range  of  cylindrical!)’ 
imploding  benchmarks  and  two  physically  plausible  test 
cases.  The  test  cases  were  primarily  treated  to  observe 
the  effects  of  a  range  cf  wire  or  gas  puff  masses.  The 
code  realistically  predicts  imploded  region  diameters, 
collapse  times,  comparative  solution  magnitudes  resulting 
from  a  variation  of  the  total  pinch  mass,  maximum  centerline 
ion  and  electron  densities,  peak  centerline  temperatures, 
and  the  radial  region  of  extent  of  radiation  emission. 

The  total  radiation  yield  is  low,  requiring  a  higher 
average  temperature  near  the  centerline  and/or  the  incor¬ 
poration  of  axially-dependent  hot  spots.  Nevertheless , 
the  presented  results  show  quantitative  realism  of  many 
parameters,  and  a  moderate  amount  of  code  alteration 
should  improve  radiation  representation. 
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Section  VIII 
FIGURE  CAP! ION'S 

Schematic  of  the  coupled  cischarge-cioae  circuit. 
KHYRAD  pressure  solution  for  the  cylindrical 
imploding  diaphragm  problem,  y  =  1.4,  1  to  1  pres¬ 
sure  ratio,  3P  finite  difference  cells. 

Lapidus  t^Ref.  II;  pressure  solution  for  the 
cylindrical  imploding  diaphragm  problem,  y  =  1.4, 

4  to  1  pressure  ratio,  199  finite  difference  cells. 
KHYRAD  density  solution,  same  conditions  as  Fig.  2. 
Lapidus  density  solution,  same  conditions  as  Fig.  3. 
KHYRAD  magnetic  field-density  ratio  solution  for 
the  imploding  cylindrical  shock  Kith  an  initially 
uniform  10“  Gauss  field  parallel  to  the  shock 
surface;  other  parameters  same  as  Fig.  2.  Infinite 
plasma  conductivity. 

KHYRAD  pressure  solution  for  the  cylindrical 
imploding  diaphragm  problem,  with  and  without 
tangential  B-field.  Infinite  plasma  conductivity 
and  other  parameters  same  as  Fig.  2. 

Simple  benchmark  timewise  voltage  distribution  for 
radiationless,  y  =  1.6"  implosions. 
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r .  Timewise  distribution  cf  radially  maximum  der.s i ty 

for  =  1.6"  Ar  with  r.o  radiation. 
riIfax (t  =  0  1  =  5  x  1018  cm  and  10  ^8  cm  ~  with 
classical,  50  times  multiplier  and  Sagdeev 
resistivities . 

10.  Radial  distribution  of  T  for  •»  =  1.6", 

e 

n^xpt  =  0  =  10" 8  cm  '  argon  at  t  =  155  nsec 

(classical  and  t  =  ITS  nsec  (50  x  multiplier  mi. 

11.  Radial  distributions  of  n,  for  the  same  conditions 
as  Fig .  10 . 


Radial  distributions 

of 

n -j  for 

Y  =  1. 

■  6"  , 

r.jax(f0i  =  5  x  101S 

cm’ 

0  Argon 

at  t 

=  ISO  nsec 

(50  x  multiplier  r) 

and 

t  =  is: 

nsec 

(Sagdeev 

15.  Radial  distribution  cf  ion  number  density,  at 

t  =  0,  measured  from  the  Physics  International 
Fi thon  Shot  No .  1440 . 

14.  Radial  distributions  of  n  and  nT  at  t  =  81.6"  and 

e  i 

101.2  nsec,  assuming  a  uniform  width,  59  cell  mesh 
with  Ar ,  T(t=0i  =  1  eY ,  the  n,  distribution  of 

J. 

Fig.  15,  and  the  voltage  distribution  of  Table  5. 

15.  Radial  distributions  of  T  for  the  conditions  of 

e 

Fig.  14. 

16.  Radial  distribution  of  radiant  power  for  the 
conditions  of  Fig.  14. 
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7 :rr.c-v, ise  di stribut icr.  of  radiant  energy  per  unit 
length  for  the  condition?  of  Fig.  14,  assuming 
both  uniform  cell  width  and  reconed  meshes. 

Mc  =  0.015  reconed  cell  widths,  at  t  =  SI. 01  and 
100.'  nsec,  for  a  39  cell  mesh  with  Ar ,  T(.t  =  0'  =  1  eV 
the  r. ,  distribution  of  Fig.  13  and  the  voltage 
distribution  of  Table  3. 

Radial  distribution  of  n  and  n,  at  t  *  SI. 01  and 

e  i 

100.'  nsec,  assuming  a  M  =  0.013  reconed  39  cell  mes 
Radial  distributions  of  T  for  the  conditions  of 
Fig.  IS. 

Radial  distribution  of  normaliced  radiant  power 
for  the  conditions  of  Fig.  IS. 

Voltage-time  profile  employed  in  the  SFLAT  and 
comparative  WHY RAD  calculations. 

Timewise  distribution  of  radiant  energy  for  6,  11, 
and  14  one  mil  Ax  wires  of  3  cm  length,  employing 
the  voltage  profile  of  Fig.  11  with  T(t*0)  =  10  eV. 
Radial  ng  and  n^  distributions  at  t  =  115. S  nsec, 
nj(t=0)  distribution  for  11  wires. 

Radial  Te  and  normaliced  radiant  power  distribu¬ 
tions  at  t  =  115. S  nsec,  T  (t*0>  distribution  for 

e 

11  wires. 

Timewise  voltage  and  current  distributions  for  a 
Maxwell  Laboratory  Blackjack  Y  pulse. 
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Racial  A R  ion  anc  electron  number  densities  at 
t  =  105."  and  118.5  nsec  for  M  =  150  ugm .  with 
the  Fig.  26  Blackjack  V  voltage  driver. 

Radial  A£  To  and  normaliced  radiant  power  distribu¬ 
tions  at  the  same  conditions  of  Fig.  2“. 

Timewise  cross  sectionally-averaged  normaliced 

radiative  power  distributions  for  M  =  50  and 

100  ugm,  with  the  Fig.  26  Blackjack  V  voltage  driver. 

Timewise  cross  sectionally-averaged  normaliced 

radiative  power  distributions  for  M  =  150  and 

500  ugm,  with  the  Fig.  26  Blackjack  V  voltage 

driver . 

Timewise  distribution  of  the  total  normaliced 
radiation  energy  for  M  =  50,  100,  150  and  300  ugm, 
using  the  Fig.  26  voltage  driver. 

Timewise  cross  sectionally-averaged  normaliced 
radiative  power  distributions  for  M  =  50  ugm  A £. 
and  Ar,  no- L  shell  line  contributions,  with  the 
Fig.  26  voltage  driver. 

Timewise  cross  sect ional lv- average  normaliced 
radiative  power  distributions  for  M  =  100  ugm  A I 
and  Ar,  no-L  shell  line  contributions,  with  the 
Fig.  26  voltage  driver. 

Timewise  distributions  of  the  total  normaliced 

radiation  energy  for  M  =  100  ugm  Ai  and  Ar,  no-L 

shell  line  contributions,  with  the  Fig.  26  voltage  driver. 


9  5 


35.  Timewise  distributions  cf  the  total  normalised 
radiation  energy  for  M  =  100  ygm  Ai ,  employing 
both  classical  r  and  50  times  anomalously  multiplie 

36.  L-shell  emission  intensity  for  M  =  150  ugm  AL  at 
t  =  105."  r.sec,  computational  radial  extent,  , 
is  0.56  cm.  0.03  £  hv  s  0.21  keV. 

3".  L-shell  emission  intensity,  same  conditions  as 

Fig.  56.  0.21s  hv  £  0.42  keV. 

38.  K-shell  emission  intensity,  same  conditions  as 

Fig.  56.  1.50  &  hv  £  2.50  k e V . 

59.  X-ray  diode  intensity  distribution,  same  conditions 

as  Fig.  36.  0.01  £  hv  £  10  keV. 

40.  L-shell  emission  intensity  for  M  =  150  ugm  At 
at  t  *  118.5  nsec,  computational  radial  extent, 

R^.,  is  0.18  cm.  0.03  s  hv  £  0.21  keV. 

41.  L-shell  emission  intensity,  same  conditions  as 

Fig.  40.  0.21  £  hv  £  0.42  keV. 

42.  K-shell  emission  intensity,  same  conditions  as 

Fig.  40.  1.50  £  hv  £  2.50  keY. 

45.  X-ray  diode  intensity  distribution,  same  conditions 

as  Fig.  40.  0.01  £  hv  £  10  keV. 

44.  U'HYRAD  program  flow  chart. 
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DIRECT  SOLUTION  OF  THE  EQUATION  OF  TRANSFER 
USING  FREQUENCY-  AND  ANGLE-AVERAGED 
PHOTON-ESCAPE  PROBABILITIES  FOR 
SPHERICAL  AND  CYLINDRICAL  GEOMETRIES 

J.  P.  Apruzese 

Science  Applications.  Inc..  McLean.  YA  Z210Z.  (J.S.A 
(Received  2  September  19R0) 

Abstract— A  previously  developed  technique  for  solvmj  the  transfer  equation  directly  by  using  frequency- 
and  angle-averaged  escape  probabilities  in  a  planar.  Doppler-broadened  medium  is  generalized  to  encom¬ 
pass  spherical  and  cylindrical  geometries,  as  well  as  a  Lorentz  opacity  profile  Two  key  elements  permit  this 
generalization  to  be  made  The  first  is  a  reciprocity  theorem  relating  the  coupling  constant  from  cell  i  to  cell 
j  to  that  from  cel!  j  to  cell  i.  The  second  is  the  use  of  a  universal  but  accurate  mean  angle  of  diffusion 

).  INTRODUCTION 

In  a  previous  paper.1  it  has  been  shown  that  line  photon-escape  probabilities  can  be  utilized  to 
obtain  exact  solutions  of  the  transfer  equation  in  a  Doppler-broadened,  plane-parallel  medium. 
The  centra;  quantity  needed  to  apply  this  technique,  the  coupling  constant  C„.  is  equal  to  the 
probability  that  a  photon  emitted  in  a  cell  of  finite  size  i  traverses  the  distance  between  cells  i 
and  j  and  is  absorbed  in  cell  j.  This  coupling  constant  is  simply  the  probability  that  a  photon 
emitted  in  cell  i  reaches  the  closest  surface  of  cell  j  minus  the  probability  that  it  reaches  the 
rear  surface  of  cell  j.  The  simplicity  of  the  plane-parallel  geometry  allows  the  frequency- 
averaged  line  escape  probability  to  be  readily  angle-averaged  as 

In  this  paper,  the  plane-parallel  formalism  is  extended  to  spherical  and  cylindrical  geometries, 
where  no  such  simple  integral  as  Eq.  (I)  applies.  In  addition,  the  technique  is  extended  to  the 
case  of  a  Lorentz  profile  by  means  of  simple,  but  highly  accurate,  analytic  fits  to  Pr  The 
technique  developed  here  may  also  be  extended  to  the  case  of  a  Voigt  profile,  provided  that 
computational!)  efficient  evaluations  of  P,  as  a  function  of  optical  depth  and  damping  constant 
are  available. 

:  CALCULATION  OF  THE  COUPLING  CONSTANTS  IN 
SPHERICAL  OR  CYLINDRICAL  GEOMETRY 

The  ease  of  determining  photon  coupling  constants  in  planar  geometry  stems  partially  from 
the  fact  that  each  cell  occupies  2 r  steradians.  as  measured  from  each  other  cell.  In  spherical  or 
c>  hndrical  geometry,  this  relation  dearly  does  not  hold  and  the  greatest  calculational  difficulties 
are  anticipated  when  computing  the  coupling  constant  from  an  outer  to  an  inner  cell,  which 
involves  complicated  integration  over  solid  angle.  Use  of  the  following  reciprocity  relation  for 
coupling  constants  overcomes  this  potential  difficulty: 

C„  =  C,jKJ  Nr  (21 

w  here  C„  and  Cr  are  the  coupling  constants  defined  in  the  introduction  and  Ar,  and  ,V,  are  equal 
to  the  total  number  of  absorbers-tatoms.  ions,  or  molecules!  in  cells  i  and  j.  respectively.  These 
coupling  constants  are  analagous  to  the  collision  probability  matrix  of  neutron  transport  theory, 
for  which  a  similar  reciprocity  relation  exists :  To  prove  Eq.  (2).  we  note  that 

K 

c„  =  2  2  c,„;n,.  <?. 

rf»  «  1  k  ■  ’ 
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where  ck„  is  the  coupling  constant  between  the  tth  atom  in  cell  j  and  the  mth  atom  in  cell  i: 
here,  the  "atom"  refers  aiso  to  ions  or  molecules.  Equation  (3 1  states  that  the  total  coupling 
constant  from  cell  j  to  cell  i  is  the  sum  of  the  coupling  constants  between  the  individual  atoms 
of  the  cells,  normalized  by  dividing  b\  the  number  of  atoms  in  the  cell  originating  the  photon. 
Since  the  coupling  constant  between  atom  k  and  atom  m  equals  that  oetween  atom  m  and  atom 
k  tsince  the  photon  path  is  the  samei  it  follows  that 


\  N, 

Q-n 


cmJK  =  2  2  cJK: 


(4) 


dividing  Eq.  (3)  by  Eq.  (4)  yields  Eq.  (2).  Equation  (2)  allows  computation  of  the  coupling 
constants  connecting  inner  to  outer  cells  only:  the  outer-to-inner  cell  coupling  constants  are 
then  obtained  from  the  reciprocity  relation. 

In  practice,  it  has  proven  possible  to  obtain  quite  accurate  results  by  tracing  the  photon- 
escape  probability  over  a  single  mean  angle  of  diffusivity  from  the  inner  to  the  outer  cell.  In 
Fig.  1.  this  process  is  graphically  illustrated  for  coupling  cell  2  to  cell  4.  The  ray  along  which  the 
P,  are  calculated  intersects  the  radius  at  the  mean  diffusivity  angle  6.  the  value  of  which  is 
discussed  below.  Since  the  escape  probabilities  must  be  averaged  over  the  originating  cell,  the 
coupling  constant  is  given  by 

C;j  =  —  [  ‘  IP, (t,  -  t) -  P,(r,  t,-  r)]dr.  (5) 

■uJo 

where  rs.  rr.  and  rf.  are  optical  depths  (see  Fig.  1)  and  P,  refers  to  the  photon-escape 
probability  along  the  chosen  mean  ray  and  is  not  angle-a\eraged.  Simple  but  quite  accurate 
expressions  for  the  single  ray  escape  probability  have  been  fitted  as  a  function  of  line  center 
optical  depth  r  for  Doppler  and  Lorentz  profiles.  For  a  Doppler  profile,  the  following 
expressions  are  used: 


P,(r)=  1/(1  -f0.65r-i-0.29r2). 

7*5.18: 

(6a) 

r  >5,18. 

(6b) 

Fig  I  The  method  used  to  obtain  coupling  constants  connecting  an  inner  cell  (2i  with  an  outer  cell  |4|  is 
illustrated  for  spherical  or  cylindrical  geometries  The  quantities  r„  ir.  and  r,  are  the  line-center  optical 
depths  along  the  indicated  segments  Escape  probabilities  are  calculated  along  the  rat  at  a  mean  angle  o 
intersecting  the  radial  at  the  center  of  the  inner  cell  The  determination  o(  f1  is  discussed  in  the  levt 
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Equation  t6a;  deviates  from  the  exact  Doppler  profile  escape  probability  in  its  domain  by  a 
maximum  of  4^  and  the  maximum  deviation  of  Eq  <6bi  is  also  4<T.  Holstein'  obtained  the 
exact  expression  corresponding  to  Eq.  (6b)  for  *  —  viz. 


which  differs  from  Eq.  (6b)  by  just  over  ?cc.  Equation  (6b)  is  utilized  rather  than  the 
mfinite-optical-depth  limit  of  Eq  (6c )  because  it  is  more  accurate  at  lower  optical  depths  where 
more  photon  coupling  occurs  and  it  joins  smoothly  at  r  =  5.lS  to  the  polynomial  expression 
given  ir.  Eq.  (6a). 

For  a  Lorentz  profile,  the  analytic  approximations  to  P,  are  also  divided  into  two  line-center 
optical  depth  domains  on  each  side  of  r  =  5.!8.  For  r«s5.18. 

P,(r)  =  l./(l.  -*•0.585727)  (7a) 

is  adopted  whereas,  for  t  >  5.18.  the  following  expression  is  used  (which  is  the  exact  limit  as 


P,(r)  =  1 'Y'(st).  (7b) 

Equations  Ca>  and  ("b)  differ  from  the  exact  escape  probability  by  a  maximum  of  897:  in  both 
cases,  the  maximum  discrepancy  occurs  in  the  vicinity  of  r  =  5. IE.  For  r<3.4  and  r  >  13.  the 
accuracy  is  better  than  39c.  The  two  expressions  join  smoothly  at  r  =  5.18. 

A  principal  question  related  to  the  use  of  this  technique  involves  the  specific  angle  P  along 
w  hich  the  rays  should  be  traced  to  obtain  the  best  coupling  constants.  These  coupling  constants 
are  directly  calculated  only  from  the  inner  cells  to  outer  cells,  as  discussed  above.  Any  photon 
emitted  from  cel!  i  whose  radius  of  curvature  is  much  less  than  that  of  cell  /  will  travel  nearly 
along  the  local  radius  w  hen  it  reaches  the  vicinity  of  the  outer  cell  j.  Therefore,  the  value  of  C,, 
obtained  is  quite  insensitive  to  the  choice  of  f  when  ceil  i  is  located  close  to  the  origin  relative 
to  cell j.  If  the  angle  fi  is  chosen  to  provide  high  accuracy  when  cells  i  and  j  have  nearly  the 
same  radius  of  curvature  (i.e..  in  the  planar  limit),  the  scheme  represented  in  Fig.  1  will 
automatically  scale  the  C„  correctly  as  cell  i  draws  closer  to  r  -  0.  Therefore,  fi  w  ill  be  chosen 
to  give  accurate  results  for  the  planar  limit,  viz. 

P,(-S)=  fp,(-)dw.  (8) 

■  cos  fi  Jq  V ' 


Since,  according  to  Eqs.  (6)  and  i").  P,ir>  depends  on  the  iine  profile,  the  best  value  of  fi  will 
also  be  profile-dependent.  In  this  paper,  the  analysis  will  be  confined  to  the  cases  of  asuophy- 
sical  and  low-to-medium  density  laboratory  plasmas  whose  line  cores  are  dominated  by  Doppler 
broadening  and  whose  line  wings  are  Loremzian.  Surprisingly  ,  the  6  which  are  chosen 
according  to  these  considerations  also  give  excellent  results  for  a  pure  Lorentz  line  (see  below ). 

The  value  of  6  which  gives  exact  results  for  large  optical  depth  is  readily  obtainable  from 
Eqs.  (6c).  (7b).  and  (8).  For  a  Doppler  profile,  letting  £  =  cos  P  where  In  r  ^  In  n- 


-\  (t  Ini  -'#2 1] 


=  f’  Mdq 
Jo  "\  ( ~  ln(  t 


r\  ( t  In  r)/,  M  2r\  l jt  In  rf^  2' 
For  a  Lorentz  profile,  we  must  have 


(9) 


(10) 
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The  error?  in  the  angie-averagcC  escape  probariit’y  a',  large  *  associated  with  a  choice  of 
£  =  cos  ••  are.  therefore,  for  Doppler  and  Lorentz  profiles 

(Doppler:  =  2,1-!  (I  la. 

and 

\P  ’ 

——  iLorentz I  =  r;\  (/II-  1.  (llbt 

M  ~ 

respectively.  Thus,  choice  of  a  £  which  is  not  exact  at  large  r  leads  to  a  constant  error  and  not 
a  constantly  increasing  error  in  P,.  Since  the  photoexcitation  of  line  radiation  is  generally 
dominated  by  absorption  of  photons  in  the  high-opacity  Doppler  core  which  are  generally 
emitted  within  a  few  optical  depths,  the  condition  of  high  accuracy  at  large  r  cannot  be  the  sole 
criterion  for  the  choice  of  £.  To  guide  this  choice,  errors  in  the  Doppler  P,  associated  with  use 
of  various  £  have  been  calculated  numerically  for  more  modest  optical  depths.  It  is  found  that, 
for  1.5  £  r<  *.  £  somewhat  larger  than  0.5  is  required  for  exact  agreement;  a  maximum  £  of 
0.54  is  required  at  -  =  3.4.  For  r<  1.5.  £  smaller  than  0.5  is  required  for  exact  agreement;  as 
r  —  0.  a!i  £  tend  toward  the  correct  P,  =  1.  In  view  of  these  results,  the  value  £  =  0.51  has  been 
adopted  for  the  following  calculations.  This  value  is  accurate  to  29f  at  -  -  *  and  is  slightly 
better  than  p.  -  0.50  for  the  important  coupling  region  1 .8  £  r  £  25.  Moreover,  even  for  a  pure 
Lorentz  profile.  £  =  0.51  is  accurate  to  better  than  l0-~  for  all  r  >  1  and.  according  to  Eq.  (1  lb), 
is  accurate  to  'Or  as  r  — No  claim  is  made  here  that  £  =0.51  is  the  best  choice  under  all 
circumstances.  It  is.  however,  a  reasonable  choice  for  our  applications  and  yields  good 
agreement  with  exact  solutions  of  several  diverse  transfer  problems.  There  is  also  a  self- 
compensating  aspect  of  the  coupling  constant  technique  for  solving  transfer  equations  which 
tends  to  reduce  the  sensitivity  of  the  results  to  the  choice  of  £.  The  source  function  in  each  cell 
is  established  by  photon  coupling  from  within  the  cell  and  from  other  cells.  If.  for  example,  the 
chosen  £  is  too  large,  the  escape  probabilities  from  other  distant  cells  to  the  cell  in  question  will 
be  overestimated,  generally  leading  to  an  overestimate  of  the  coupling.  However,  the  overes¬ 
timate  of  P,  from  the  local  cell  partially  compensates  for  this  effect  by  reducing  the  local 
coupling  which  is  proportional  to  (I  -  P,\.  As  is  shown  in  Ref.  1.  small  errors  in  the  coupling 
matrix  do  not  result  in  large  source-function  errors. 

As  a  precaution  and  to  ensure  consistency  ,  the  X  coupling  constants  from  cell  /  to  other 
cells  j  are  renormalized,  as  will  now  be  described.  First,  the  probability  of  escape  from  the 
entire  plasma  from  cell  i  is  calculated  along  the  chosen  6.  This  quantity.  P,.  must  be  equal  to 

S 

1  -  S  C„.  It  can  be  found  to  be  slightly  different  because  of  numerical  noise:  therefore  the  C,  are 
multiplied  by  the  factor  needed  to  make  this  sum  over  j  equal  to  (1  -  P,  i. 

3  COMPARISON  OF  RESULTS  WITH  EXACT  SOLUTIONS 
The  most  meaningful  test  of  our  formalism  is  the  following:  how  closely  are  we  able  to 
reproduce  exact  solutions  of  basic  transfer  problems  in  spherical  or  cylindrical  geometry?  The 
technique  of  using  the  coupling  constants,  once  obtained,  to  solve  for  the  steady-state  line 
source  function  is  described  in  Ref.  I.  Briefly,  the  equation  for  the  steady -state  upper  level 
population  of  a  two-level  atom  in  cell  i  is  written  as 

=  o  =  \,,H ,  -  X ‘V.A.-C,  -  N.,1 .4..  -  D, I-  (12. 

where  .\„  and  X,,  are  the  upper  and  lower  level  populations  of  cell  i.  respectively.  A,  is  the 
spontaneous  transition  probability  (sec'1).  C,.  is  the  coupling  constant  from  cell  i  to  cell  i  as 
previously  defined,  and  W,  and  D,  are  the  total  collistonal  population  and  depopulation  rates, 
respectively,  of  the  upper  level  in  cell  i.  For  the  cases  studied,  the  absorption  coefficient  and 
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hcnvc  the  C  .-.re  fixed  A  fixed  absorption  coefficient  also  allows  the  lower  lex  el  popu.utions  .Y 
to  re  written  in  term'-  of  ihe  upper  icxei  population  \„  for  a  two-ievel  atom  xir.ee 

k  r,  S,B„  -  S.B. . 

wnerc  ihe  B  are  the  stimulated  absorption  and  emission  probabilities.  Solving  Eq.  « 12 >  for  the 
two-ies el-atom,  thus  reduces  to  an  mxersior.  of  an  V  *  X  matrix.  The  requirements  imposed  on 
the  W.  D  and  .4.  stem  from  the  fact  that  the  comparison  cases  specify  the  Pianck  function 
spatially  and  aiso  the  quenching  parameter  Pt  w  hich  is  the  ratio  of  collisional  to  total  decay  of 
the  upper  state.  The  W  .  D  and  A„.  of  Eq  <  12 1  are  chosen  to  be  compatible  w  ith  the  specified 
Planck  function  and  quenching  parameter  The  C  are  calculable  from  the  fixed  absorption 
coefficient. 

The  results  of  comparisons  with  previously  published  solutions  to  basic  transfer  problems 
are  presented  in  Figs.  2-4.  In  Fig.  2.  comparison  is  made  with  the  hollow -sphere  calculations  of 
Kunasz  and  Hummer.4  The  cylindrical  calculations  of  Avery  el  al:  are  repeated  with  the 
present  method  and  the  results  exhibited  in  Fig.  3.  The  comparisons  of  Fig  3  are  confined  to 
cases  where  the  cylinders  are  long  compared  to  their  widths,  since  the  present  techniques  apply 
to  infinite  cylinders.  Finally,  since  the  mean  diffusion  angle  used  in  these  calculations  is  most 
accurate  for  a  Doppler  profile,  a  comparison  with  pure  Lorentz  profile  solutions  obtained  by 
Avretl  and  Hummer'  is  presented  in  Fig.  4  to  determine  the  applicability  of  p.  =  0.51  to  Lorentz 
and  hence  to  Voigt  profiles. 

It  is  self-evident  from  the  three  figures  that  the  method  yields  very  good  accuracy.  The 
maximum  deviation  of  the  calculated  source  function  from  that  given  by  accurate  solutions  is 
259  for  the  cylindrical  and  spherical  cases,  with  most  of  the  curves,  lying  within  159c  of  the 
previously  obtained  values.  For  the  planar  Lorentz  profile,  ihe  maximum  error  is  109.  Although 
exact  solutions  are  clearly  not  obtained  by  using  this  formalism,  the  accuracy  is  more  than 
acceptable  considering  that  the  technique  obviates  entirely  the  need  for  either  frequency-  or 
angle-integration.  As  might  be  expected,  the  largest  number  of  cells  required  to  obtain  this 
accuracy  (75 1  was  needed  for  the  stringent  lest  of  the  Kunasz-Hummer  cases,  where  both 
the  absorption  coefficient  and  the  Planck  function  vary  as  r":.  The  smallest  number  of  cells  was 


Fig  2  Sieadv -stale  source  function*  obtained  n>  use  of  ihe  presently  desenped  method  idoiicd  line )  are 
.orr.pLf-.-c  with  ev^ei  -phencal  solutions  obtained  in  Ref  A  h\  Kunast*  ,ind  Humme'  noiid  line •  in  obtaining 
the  prefer.*  jprroMnjti  solutions  we  used  *5  celis  ioganthmicalh  spaced  in  -  Tht  radiu*  of  the  inner 
boundu"*  of  each-  sphere  i«  J  and  each  cur>e  i*  labeled  hi  ihe  radius  of  the  ouier  rvundan  The 

fun*!*'  r  ,:nc  j^orpiion  coefficient  vary  a*  indicated  A  pun-  Doppier  profit-  is  assumed  with  (me 
center  opt::;.:  depth  W'\  r  =  5 l«i.Z 
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Fig  3  Steady -stale  source  (unctions  obtained  by  the  present  method  (dotted  lines )  2re  compared  with  the 
Monte  Carlo  cylindrical  solutions  obtained  in  Ref  5  isolid  iinest  The  curses  ate  differentiated  by  line 
center  optical  depth  along  the  cylinder  radius.  Ir.  obtaining  the  present  approximate  solutions.  30  cells 
logarithmically  spaced  in  r  were  used,  with  selected  cells  plotted  The  iines  connecting  the  points  are 
intended  to  guide  the  eye  rather  than  to  suggest  interpolated  saiues 
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Fig  4.  Steady-state  source  (unctions  obtained  by  the  present  method  (circ)cs i  are  compared  with  the  e\act 
planar  Lorentz-profile  solutions  obtained  in  Ref  b  (solid  linesl  The  Planck  (unction  is  assumed  to  be 
spatially  uniform  and  equal  to  I.  The  number  of  cells  (ioganthmically  spaced  in  ri  presently  used  is  10.  10. 
35.  and  40.  ir  ascending  order  of  optical  depth  The  value  of  the  source  function  at  the  center  of  each  cel! 
is  indicated  bs  the  circles,  except  for  r  =  1.593  *  10'. for  which  every  other  celi  plus  the  innermost  cell  are  plotted 


used  (10)  in  calculating  the  low-optical-deplh  planar  Lorentz  cases  for  which  the  Planck 
function  was  spatially  uniform. 

4  SUMMARY  AND  CONCLUSIONS 

It  has  been  shown  that  the  use  of  escape  probability  techniques  to  ootain  solution1  to  the 
radiative  transfer  equation  ma>  be  extended  from  planar  to  cylindrical  and  spherical 
geometries  The  formalism  for  line  transfer  requires  that  the  line-profile-averaged  photon 
escape  probability  along  a  ra\  be  known  as  a  function  of  optical  depth  Convenient  ar.aotic  fits 


Direct  soiuiK'r.  of  the  equation  of  ttansfet 


to  these  quaniuies  are  given  (or  the  cases,  of  pure  Doppler  and  pure  Dorentz  line  profiles.  The 
technique  presented  is  quite  advantageous  numerically  in  that  it  eliminates  frequency  and  angle 
integrations  at  small  cost  in  accuracy.  It  is  also  readily  adaptable  to  different  physical  situations. 
For  instance,  for  high-density  plasmas  for  which  Stark  broadening  dominates  the  line  profile,  the 
escape  probability  formulae  given  by  Weishei’.'  could  be  readily  employed 

Perhaps  the  mam  drawback  to  our  approach  (and  the  principal  area  for  further  work  on  this 
technique1  is  that  profile-averaged  escape  probabilities  are  most  easily  obtained  and  used  for 
spatially  constant  line  profiles.  When  density  and/or  temperature  gradients  are  present  in  such 
degree  as  to  affect  P,(rl  significantly,  the  accuracy  of  this  approach  will  deteriorate.  Obtaining 
tits  to  P,i- )  in  the  presence  of  such  gradients  is  the  only  way  to  restore  accuracy. 

Acknoviedpmtnt— This  work  was  supported  by  the  Naval  Research  Laboratory 
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Bv  self-consistently  calculating  x-rav  spectra  from  first  principles,  we  have  delineated  the  relationships  between  the 
spectrum  and  the  state  of  compression  and  heating  of  a  neon  plasma  in  detail  A  coUtsionai-radiame  mode1 
including  Stark  line  profiles  is  used  to  determine  the  highly  ionized  high-density  states  of  neon  One  of  our 
calculated  spectra  reproduces  remarkably  well  an  experimental  spectrum  obtained  from  laser  implosion  by  Vaakob: 
cl  a!,  and  indicates  compression  conditions  significantly  different  from  those  obtained  assuming  the  validity  of  locai 
thermody  namic  equilibrium  Implications  of  our  calculations  for  spectroscopic  diagnoses  of  fusion  plasma,  .u 
discussed 


1  INTRODUCTION 

One  of  the  most  widely  used  techniques  to  diag¬ 
nose  the  properties  of  transient  fusionlike  plasmas 
is  to  analyze  their  x-ray  emissions.  In  an  example 
of  the  use  of  such  methods.  Yaakobi  ct  nl have 
measured  the  A' -series  emission  spectrum  from  a 
neon  plasma  produced  by  laser  implosion  of  a 
glass  microballoon  filled  with  neon  at  8.6  atm 
pressure.  Using  a  methodology  that  assumed  the 
ionization  conditions  of  the  plasma  and  that  ac¬ 
counted  for  opacity  and  Stark  broadening  of  the 
Lyman-o  and  Stark  broadening  of  the  Lyman-.! 
and  Lyman-}  lines,  they  ascribed  a  density  of 
0.26  g  cnr  (ion  density.  7.8  >  10!1  cm*si  to  the 
neon  emission  region.  From  the  measured  line 
intensity  ratios  /(Ly-i)  7(ls:-ls3/>)  and/(Lyoi 
1{ Lyj),  an  electron  temperature  of  300  eV  was 
inferred  in  the  context  of  an  LTE  plasma  model 
whose  line  radiation  was  assumed  to  be  a  Planck 
function — saturated.  In  this  paper,  we  show  that  a 
theoretically  sell-consistent,  steady-state  analysis 
of  all  the  main  features  of  such  spectra  can  be 
carried  out  that  is  based  on  a  first -principles  cal¬ 
culation  in  which  the  only  free  parameters  are  the 
average  plasma  properties  of  temperature,  densi¬ 
ty.  size,  and  velocity  profile.  This  analysis  leads 
to  a  set  of  general  principles  governing  the  rela¬ 
tionships  of  high  density,  optically  thick  plasma 
properties  to  their  emitted  spectra.  As  a  specific 
example  of  the  application  of  such  principles,  we 
demonstrate  that  our  calculations  reproduce  well 
the  neon  experimental  spectrum  discussed  above 
and  lead  to  a  diagnosis  of  significantly  lower  den- 
sitv  and  higher  temperature.  Furthermore,  our 
investigation  allows  an  exposition  of  how  and  why 
a  self-consistent  physical  model  reproduces  an 
experimental  spectrum,  while  a  simpler  bui  in¬ 
appropriate  model  such  as  LTE  may  lead  to  mis¬ 
calculation  and  or  misinterpretation  of  the  spec¬ 
trum. 


11  PLASMA  MODEL 

For  spectrum  calculations  we  have  employed  a 
detailed  multistate,  multilevel  model  of  ionized 
neon  in  the  context  of  coilisional-radiative  equilib¬ 
rium  and  an  assumed  spherical  geometry.  Since 
virtually  all  of  the  neon  at  peak  compression  con¬ 
sists  of  lithiumlike  or  more  highly  stripped  ions, 
the  model  includes  only  the  ground  states  for 
Ne  I-Ne  VII.  For  lithiumlike  Ne  vm.  the  atomic 
model  includes  the  excited  states  l-:-2/>.  -3-.  -3/>. 
-3 and  -4d.  For  heliumlike  NeIX.  we  have  the 
excited-state  manifold  ls2s 's.  l>2/>  'I1 .  1.-2  p'P. 
and  the  u  =  3  and  u  =  I  singlet  states.  For  Ne  X. 
n  ~2 .  3 .  4 .  and  5  as  well  as  n  =  1  are  included. 

All  ol  the  iines  appearing  m  the  Rochester 
spectrum  tLyo.  p.  and}  .  NeIX  l.«2-ls 3/>  'P. 
U:-ls4f>'F.  as  well  as  ls!-l-2 f>'P)  are  calcu¬ 
lated  by  sell-consistently  solving  fully  coupled 
nonlinear  radiation  transport  and  rate  equations." 
Most  importantly,  the  calculation  also  includes 
the  sell-consistently  computed  effects  of  Doppler. 
Stark/  and  mass-motional  line  broadening.  Addi¬ 
tionally.  the  calculated  spectra  are  convolved  with 
1.5  eV  of  experimental  (Gaussian)  broadening/ 
However,  m  most  cases,  as  pointed  out  by  Yaakobi 
rl  nl ..  the  experimental  broadening  has  little  or  no 
effect  due  to  the  large  Stark  linewidths.  The  rale 
equations  and  the  radiative  transport  equation  are 
also  used  to  solve  in  detail  for  the  continuum  radi¬ 
ation  arising  from  recombination  from  the  bare 
nucleus  to  the  n  -  1  and  2  states  of  Ne  X.  and  from 
Ne  X  n  =  1  to  Ne  IX  1-"  *S  and  1-2 p  3/V 

With  the  exception  of  electron  collisional  excita¬ 
tion  rate  coefficients  and  spontaneous  decay  rates, 
the  methods  ol  calculating  the  various  rate  coeffi¬ 
cients  used  in  this  study  have  been  described  in 
detail  in  previous  papers. i‘  and  only  brie!  telet- 
ence  w  ill  be  made  lo  these.  Everv  slate  is  cou¬ 
pled  to  the  next  ienergeticallv>  highest  urnuml 
slate,  by  collisional  ionization  .  1  collis  i  re- 
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combination.  Additionally,  pbotoionization  and 
radiative  recombination  are  computed  in  detail  lor 
tne  states  mentioned  above.  The  collisionai  ion¬ 
ization  rates  are  calculated  by  Seaton's  prescrip¬ 
tion,'  and  the  p  lotoionization  cross  sections  are 
calculated  in  the  hvdrogenic  approximation7  using 
effective  free-bound  Gaunt  factors.8  Collisional- 
and  total-radiative  recombination  rates  are  then 
calculated  as  the  detailed  balance  of  these.  How¬ 
ever,  the  details  of  the  frequency  dependence  of 
the  hydrogenic  photoionization  cross  sections  are 
replaced  by  a  simpler  exponential  profile  for  the 
recombination  process.5  In  addition,  adjacent 
ground  stales  are  coupled  by  dielectronic  recombi¬ 
nation8  found  by  summing  capture  rates  over  a 
manifold  of  states  above  the  ionization  limit  which 
then  decay  via  cascade  to  the  ground  state. 

Excited  levels  of  a  given  ion  are  coupled  to  the 
other  excited  states  in  the  model  and  to  the  ground 
state  by  electron  collisionai  excitation  and  de- 
excitation.  and  by  spontaneous  radiative  decay. 
Collisionai  couplings  include  forbidden  and  spin- 
flip  transitions  as  well  as  those  which  are  dipole 
allowed.  ForNeX,  the  Coulomb-Born  approxima¬ 
tion17  was  used,  while  a  distorted-wave  calculation 
with  exchange"  was  used  to  calculate  the  coeffi¬ 
cients  for  Nevm  andNeiX.  Comparisons  of  these 
two  methods  have  been  made  for  several  transi¬ 
tions  in  hydrogenlike  A1  XIII ,!!  and  excellent  agree¬ 
ment  was  obtained. 

As  mentioned  above,  line  photon  reabsorption  is 
taken  into  account  by  solving  the  radiation  trans¬ 
port  equation  for  six  selected  resonance  lines, 
five  of  which  appear  in  the  published  Rochester 
spectrum1  (taken  by  Yaakobi  el  al.)  and  are  gener¬ 
ally  used  as  diagnostics  for  neon-  or  argon-filled 
pellet  targets.  The  present  model  neglects  the 
effect  of  radiation  reabsorption  for  lines  which 
couple  two  excited  states,  even  though  such  lines 
are  optically  thick  in  some  of  the  cases  discussed 
beiow.  This  simplification,  however,  results  in 
no  significant  error  in  these  cases,  where  the 
density  is  so  high  that  collisionai  rates  coupling 
excited  states  vastly  exceed  radiative  rates  and 
thus  dominate  the  cross  coupling  of  these  levels. 
Unless  these  lines  appear  in  the  spectrum  to  be 
analyzed,  it  is  not  necessary  to  perform  the  radi¬ 
ation  transport  calculations  for  them.  The  same 
'principle  does  not  apply  for  the  resonance  lines 
dominating  most  spectra,  where  the  reduced  col¬ 
lision  rales  and  higher  radiative  rates  result  in 
comparable  importance  for  collisions  and  radia¬ 
tion.  or  even  radiation  dominance  in  the  couplings. 

Details  of  the  radiation  transport  and  photocou¬ 
pling  calculations  are  given  elsewhere5  ,IS  and  will 
not  be  elaborated  here.  However,  ar.  additional 
feature— the  calculation  of  the  effects  of  the  dif- 
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ferent  Doppler  shifts  caused  by  tne  inhomogeneous 
velocity  profile  of  a  plasma  implosion  on  tne 
emitted  line  profiles — has  been  included  in  tne 
radiation  transport  algorithm.  Tnese  efiects  will 
be  fully  discussed  in  an  upcoming  paper,  but  in 
the  particular  results  of  this  paper  it  served  us 
mainly  to  confirm  what  nad  been  suspected,  that 
even  at  implosion  velocities  of  2  MO1  cm  sec", u 
the  effects  are  small  due  to  the  large  Stark  widths 
of  the  lines  at  high  density.  For  lower  density 
plasmas  the  effects  can  be  quite  significant  but  it 
is  of  little  concern  here. 

In  the  analysis  detailed  below,  spectra  calculated 
for  steady-state  plasma  conditions  are  used  to  in¬ 
terpret  an  experimental,  time -integrated  spec¬ 
trum.  While  this  approach  may  at  first  seem 
questionable,  it  was  justified  for  the  cases  con¬ 
sidered  here — short-pulse  Nd  laser  drivers — by 
Yaakobi  et  al.15  in  their  letter  on  argon-filled  glass 
microballoon  experiments.  They  found  from  ex¬ 
tensive  hydrodynamical  simulations  of  spherical 
target  compressions  that  the  most  energetic  emis¬ 
sion  of  hydrogenlike  lines  occurs  during  a  10-psec 
time  interval  during  which  the  plasma  properties 
change  little,  spatially  or  temporally.  These  theo¬ 
retical  results,  coupled  with  the  fact  that  our 
analysis  obtains  striking  agreement  with  experi¬ 
ment,  indicate  that  this  assumption  can  be  pro¬ 
ductively  employed  under  these  conditions  to  inter¬ 
pret  the  emission  properties  of  plasma  cores  under 
spherical  compression. 

III.  COMPRESSION'  DIAGNOSTICS 
WITH  CALCULATED  SPECTRA 

It  is  widely  recognized  that  pellet  compression 
along  a  low  adiabat  is  most  desirable  in  achieving 
maximum  density  prior  to  thermonuclear  ignition. 
Hence,  a  major  goal  of  laser-driven  pellet  com¬ 
pression  experiments  is  to  obtain  the  highest  pos¬ 
sible  implosive  drive  while  minimizing  preheat. 
Knowledge  of  the  degree  to  which  this  goal  is  met 
experimentally  may  be  obtained  from  a  variety  of 
techniques,  and  one  of  the  prime  tools  available 
is  spectroscopy.  For  this  purpose,  one  must  know 
the  sensitivity  of  the  emitted  spectrum  to  the 
plasma's  density  and  temperature.  Figure  1  pre¬ 
sents  three  neon  spectra  calculated  for  a  spherical 
plasma  in  CRE.  They  illustrate  the  evolution  of 
the  emitted  spectrum  as  the  same  mass  of  neon  is 
progressively  compressed  to  higher  densities  and 
low-er  temperatures  (i.e..  along  lower  adiabatsi. 
The  lower  the  adiabat.  the  closer  to  tne  lop  of  the 
figure  the  spectrum  appears.  Increased  density  is 
primarily  reflected  in  the  widening  of  the  lines  due 
both  to  larger  intrinsic  Stark  widths  and  possible 
opacity  broadening.  This  effect  has  been  utilized 
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FIG.  1.  Three  calculated  spectra  are  shown  which  il¬ 
lustrate.  Irom  bottom  to  top,  the  changes  in  spectral 
signature  which  occur  as  a  constant  neon  mass  is  com¬ 
pressed  to  higher  densities  and  colder  temperatures. 

The  spectra  are  calculated  assuming  that  the  neon 
plasma  is  in  collisional-radiative  equilibrium  (CREi. 

The  exterior  pellet  temperatures  decrease  from  those  of 
the  core  in  accordance  with  the  gradient  descrilied  in 
the  text.  Oily  the  middle  spectrum  includes  the  (negli¬ 
gible!  pumping  effect  of  the  Na  v  impurity  line  at  1.13 
keV. 


for  argon  compression  diagnostics  in  the  case  if 
optically  thin  nonopacity  broadened  Stark  lines." 
The  line  profiles  las  opposed  to  line  intensities) 
are  relatively  insensitive  to  temperature.1'  How¬ 
ever.  effects  of  opacity  broadening  are  specifically 
illustrated  below  in  the  interpretation  of  the  Roch¬ 
ester  neon  experimental  spectrum1  taken  by  Yaakobi 
r  /  al . 

The  pronounced  temperature  sensitivity  of  the 
emitted  spectra  is  also  apparent  in  Fig.  1.  As  the 
temperature  is  decreased  from  600  to  450  to  250 
eV .  the  intensity  of  the  helium-like  ls:-lx3 p  '/' 
line  at  1.074  keV  increases  dramatically  relative 
to  the  hydrogen-like  NeX  lines.  In  the  Rochester 
spectrum,  the  NelX  l.s- - Is 4/>  'P  line  (at  1.127  keV) 
was  obscured  by  the  impurity  line  Nax  ls:-l.s2/>  '/J. 
which  is  equal  in  wavelength  to  1  part  in  10s.  In 
fact  this  sodium  line  will  optically  pump  the  NelX 
level  populations  to  some  degree  due  to  this  close 
resonance.  In  the  middle  spectrum  of  Fig.  1  the 
effect  of  this  impurity  line  has  been  included — the 
intensity  having  been  derived  from  our  computed 
fit  to  the  Rochester  spectrum'  (taken  by  Yaakobi 
i  /  al .)  in  Fig.  2.  We  included  the  sodium  pumping 
line  by  assuming  that  the  flux  directed  inward 
from  the  nucroballoon  toward  the  neon  filler  is 
the  same  as  the  inferred  flux  directed  outward. 
However,  by  calculating  the  emitted  neon  spectrum 
without  the  pumping  influence  of  the  impurity  line, 
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FIG.  2.  The  fop  spectrum  was  obtained  by  Yaatiobi  rt 
at.  at  the  University  of  Rochester  from  a  laser-imploded 
microballoon  filled  with  s.O  atmospheres  ol  neon.  The 
bottom  two  spectra  were  calculated  assuming  CRK.  The 
middle  spec'rum  provides  the  >est  ill  to  the  top  one  and 
arises  irom  a  plasma  of  T,.=  3x,V  eV,  1  =  11.4  urn  at 

'  J  (»-•  cm''1.  The  bottom  spectrum  would  lie  emitted 
Dv  a  plasma  of  T.«300  eV.  r»  pm,  al  \t-  - 1 l" 

cm'3. 


we  have  determined  that  the  effeci  is  negligible  for 
these  dense  plasmas  with  high  collisional  cross¬ 
coupling  rales  between  excited  stales.  Therefore, 
the  line  has  been  omitted  in  some  ol  the  calculated 
spectra  presemed  in  Figs.  1-3.  since  11s  effeci  on 
the  level  populations  is  negligible. 

Equally  apparent  as  the  helium  line  intensity 
change  is  the  relative  increase  in  the  continuum 
radiation  arising  from  recombination  10  the  ground 
slate  of  Ne  IX  as  the  temperature  decreases.  In¬ 
deed.  at  a  temperature  of  250  eV  the  Lyg  and  Ly, 
lines  which  sit  on  fop  of  this  continuum  feature 
fwhose  edge  lies  at  1.196  keV)  are  more  pruperlv 
regarded  as  “  line-continuum  lealures"  rather  than 
pure  lines.  The  appearance  ol  tins  strung  continu¬ 
um  edge  can  place  strong  constraints  on  the  diag¬ 
nosed  temperature.  Also,  the  continuum  opacity 
may  significantly  aflect  both  the  shape  anti  strength 
of  the  observed  features  arising  from  bound  level 
transitions.  Similarly,  ionization  bv  line  pho¬ 
tons — especially  important  when  the  lines  are  very 
wide  due  to  Siark  effects  — ma\  sigmlicanth  aflect 
(he  ionization  balance  of  the  plasma. 

Having  noted  Hie  spectral  differences  arising 
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FIG.  3.  The  three  calculated  spectra  shown  represent 
the  results  of  employing  commonly  used,  but  inaccurate, 
approximations.  The  top  spectrum  would  be  emitted  by 
a  neon  plasma  with  pure  Voigt  line  profiles  whose  condi¬ 
tions  are  those  of  the  "best  fit"  of  Fig.  2.  The  middle 
spectrum  maintains  the  same  physical  conditions  as  the 
top,  but  Stark  line  profiles  and  LTE  are  assumed.  The 
bottom  spectrum  arises  from  the  same  mass  of  plasma 
assumed  to  be  in  LTE  but  with  a  temperature  of  300  eV 
and  a  radius  of  10  pm.  The  sodium  imparity  line  at 
1.13  seV  is  included  in  these  spectra. 


from  changes  in  density  and  temperature,  we  in¬ 
vestigated  the  possibility  of  obtaining  a  fit  to  the 
Rochester  experimental  spectrum1  using  our  first- 
principles  ionization  calculation,  in  which  only 
temperature,  density,  size,  and  velocity  variations 
in  plasma  conditions  could  be  made. 

The  results  of  this  effort  are  graphically  sum¬ 
marized  in  Fig.  2,  where  the  Rochester  experi¬ 
mental  spectrum  (taken  from  Ref.  1)  and  two  cal¬ 
culated  spectra  are  displayed.  The  middle  spec¬ 
trum  of  Fig.  2  represents  an  excellent  fit  to  the 
experimental  spectrum  displayed  at  the  top  of  the 
figure.  This  "  best-fit"  calculated  spectrum  would 
arise  from  a  spherical  neon  plasma  of  ion  density 
4.5xio!t  cm':  and  radius  11.4  pm,  having  an  elec¬ 
tron  temperature  of  385  eV  from  r  =  0  to  r=8.8 
pm  that  then  decreased  smoothly  to  296  eV  at  the 
surface.  This  gradient  was  inferred  in  a  previous 
analysis  of  the  degree  of  line  self-reversals  in 
such  imploded  plasmas.11  and  our  present  results 
support  this  gradient  inference.  Also  shown  in 
Fig.  2  i at  the  bottom)  is  the  spectrum  which  would 
be  emitted  by  a  neon  plasma  characterized  by  the 
conditions  originally  inferred  from  the  LTE  analy¬ 
sis  in  Ref.  1.  These  conditions  are  T,  =  300  eV 
(constar.il.  .V,  =  7.B'  10:i  cm",  and  v  =  9.5  pm. 

The  need  for  a  higher  temperature  and  lower 
density  than  previously  diagnosed  is  most  readily- 


understood  by  comparing  the  predicted  300-eV 
spectrum  to  the  experimental  one.  The  presence 
of  the  prominent  continuum  feature — the  Ne  IX 
ground-state  ionization  edge— in  the  300-eV  pre¬ 
dicted  spectrum  is  not  matched  by  a  similar  fea¬ 
ture  in  the  experimental  one.  An  explanation  for 
its  absence  is  that  the  experimental  plasma  con¬ 
tains  less  NeK  and  is  therefore  at  a  higher  tem¬ 
perature,  where  more  of  the  heliumlikc  neon  is 
ionized.  This  degree  of  reduction  of  the  continuum 
feature  as  the  temperature  increases  is  clearly 
evident  in  the  theoretical  spectra  of  Fig.  1,  which 
served  as  a  guide  in  obtaining  the  Fig.  2  fit  to  the 
Rochester  experiment.  A  core  temperature  of  385 
eV  was  needed  to  lower  the  calculated  feature 
which  is  weak  or  nonexistent  in  the  experimental 
spectrum.  Moreover,  temperatures  higher  than 
385  eV  were  found  to  produce  a  Neix  lss-ls3p‘P 
line  which  was  too  weak  to  match  the  experimental 
line.  Even  at  7  =  400  eV,  the  Fig.  2  comparison  of 
calculated  and  experimental  spectra  deteriorated 
noticeably  in  this  regard.  Any  temperature  below 
385  eV  raises  the  helium  continuum  again  produc¬ 
ing  a  significant  deterioration  in  the  agreement 
even  at  our  diagnosed  density  of  4.5  xlO21.  Cur 
ability  to  zero  in  on  385  eV  as  the  likely  tem¬ 
perature  in  this  manner  strikingly  illustrates  both 
the  need  and  the  value  of  utilizing  all  of  the  infor¬ 
mation  in  the  spectrum. 

The  Rochester  analysis  yielding  Tc  -  300  eV  was 
based  upon  Planckian  line  ratios  generated  by  an 
LTE  plasma.  By  contrast  we  find  that  a  neon  plas¬ 
ma  of  7^  =  300  eV,  .Y(.  =  7.8  xlO21  cm'5,  and  r  =  9.5 
pm  has  not  reached  LTE.  Of  the  ionic  species, 
only  the  NeX  concentration  is  close  to  its  LTE 
value.  NeXJ  (bare  nucleus)  is  about  half  of  what 
would  be  expected  in  LTE,  while  heliumlike  Neix 
has  twice  its  LTE  concentration.  Mo-eover,  the 
excited  states  of  NeX,  whose  radiative  decay  dom¬ 
inates  the  spectrum,  are  present  to  a  far  smaller 
degree  than  would  be  the  case  in  LTE.  The  LTE 
excited-to- ground- state  ratios  of  NeX  n  =  2.  3.  4. 
and  5  are,  respectively,  0.13,  0.16,  0.23.  and 
0.32.  We  find  that  the  highest  ratios  occur  in  the 
300-eV  model  plasma  at  >-  =  0.  where  they  are 
0.05,  0.06,  0.09,  and  0.13,  respectively.  In  effect, 
the  ionization-excitation  state  of  the  plasma  lags 
behind  what  is  required  for  LTE.  At  an  ion  density 
of  7.8  xlO21  cm'3,  collisional  processes  dominate 
the  coupling  of  the  NeX  excited  states  to  each 
other.  Hence,  these  ratios  of  populations  of  the 
states  lying  above  >i  =  2  to  that  of  u  =  2  do  corres¬ 
pond  to  the  required  LTE  ratios.  The  NeX  excited 
states  are  also  in  LTE  with  respect  to  the  ground 
state  lying  above  (the  bare  nucleus)  due  t  >  the 
strength  of  collisional  ionization  and  remmbination 
processes.  However,  collisional  excitation  from 
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and  de-excitatiup.  to  the  NeX  ground  state  from  the 
excited  states  is  generally  at  least  two  orders  of 
magnitude  weaker  than  collisional  cross  couplings 
among  the  excited  states.  Collisional  deexcita¬ 
tion  from  the  excited  states  to  the  ground  state  is 
also  only  about  one-third  as  likely  as  spontaneous 
radiative  decay  to  the  ground  state.  Collisional 
processes  are  therefore  not  strong  enough  to  en¬ 
force  LTE  between  the  ground  and  excited  states 
of  the  sai,  e  species— NeX.  Moreover,  due  to  the 
large  Stark  widths,  the  optical  depths  of  Lyp.  >, 
etc.,  are  2.2  or  lower,  as  described  below.  Since 
these  modest  opacities  do  not  permit  multiple 
photon  scattering,  the  states  h=3  and  above  are 
weakly  photoexcited  and  remain  far  from  LTE 
with  respect  to  n  =  1.  However,  since  the  optical 
depth  of  Lya  is  -200,  it  might  at  first  be  sus¬ 
pected  that  at  least  n  =  2  would  be  in  LTE  with  re¬ 
spect  to  the  ground  state.  The  fact  that  this  is  not 
the  case  may  be  understood  by  applying  the  condi¬ 
tions  required  for  LTE  to  the  rate  equations  af¬ 
fecting  the  two  states  n  =  l  and  2.  For  simplicity 
we  first  consider  the  two- level-atom  case  in  an 
optically  thick  plasma  characterized  by  a  line 
photon  escape  probability  P,.  We  let  H2J  stand  for 
the  collisional  deexcitation  rate  coefficient  and 
Ajj  stand  for  the  spontaneous  emission  coefficient. 
The  steady-state  condition  requires 

ill 

where  photon  reabsorption  is  accounted  for  by 
diluting  the  spontaneous  emission  coefficient  by  a 
factor  equal  to  the  mean  photon  escape  probability. 
From  Eq.  ( 11, 


ly  thick  interior  based  on  the  abuve  two- level- 
atom  analysis.  However,  there  are  three  separate 
effects  which  in  this  case  prevent  LTE  conditions 
from  being  established  even  in  the  thick  plasma 
interior.  First,  even  though  the  optical  depth  of 
Lya  is  -200,  the  line  wings  are  so  prominent  at 
this  high  density  that  the  mean  escape  probability 
is  found  numerically  to  be  0.0b.  Thus.  A2.P 
~0.15  H21  in  the  interior.  According  to  Eq.  2i 
this  alone  will  reduce  the  ,V2  .V,  ratio  to  0.67 
times  that  of  LTE.  The  other  twu  eflects  which 
bring  the  actual  ratio  to  only  0.3B  times  the  LTE 
ratio  are,  first,  a  radiative  leak  in  the  Lyp  line, 
where  ts2,  and,  second,  another  photon  escape 
mechanism  in  the  optically  grey  recombination 
continuum  from  the  bare  nucleus  to  the  ground 
state  of  NeX.  A  multilevel  analysis  is  necessary 
to  clarify  the  detailed  contributions  of  the  latter 
two  mechanisms. 

An  accurate  analytic  treatment  is  possible  in  this 
case  because  the  very  large  collisional  rates  coup¬ 
ling  the  excited  states  (n  =  2,  3,  4,  5,  and  NeXl) 
force  the  ratios  of  these  states  to  each  other  to 
equal  those  of  LTE,  as  confirmed  by  the  calcu¬ 
lated  population  densities.  Let  £  denote  the  LTE 
ratio  of  state  i  to  state  i.  Then,  the  equation  for 
the  steady- state  coupling  of  NeX  u  =  1  to  the  over- 
lying  states  becomes 

A'.(E  + 

o' 

I*  3 

+  M«WVv^l)-  4) 


V  I  .  ,r 

N,  W2l> 


(2) 


For  .V2  .V.  to  have  the  LTE  population  ratio  of 
these  two  levels,  ( ,'2  we  must  have 

A2lP,.<  U21  .  (3) 

Therefore,  even  a  very  optically  thick  plasma 
can  be  far  out  of  LTE  if  A21  U'21  to  the  extent 

necessary  to  offset  a  small  P  .  That  is,  the  mere 
existence  of  a  small  P  due  to  large  optical  depth 
does  not  ensure  the  validity  of  Eq.  13).  Therefore, 
the  lact  that  Lya  is  quite  thick  in  this  case  does 
not  in  itself  justify  an  LTE  assumption.  For  the 
plasma  conditions  under  consideration  here. 

A2J  s  2.5H21.  Therefore.  LTE  will  not  be  valid 
near  the  surface  of  the  plasma,  where  P  *0.5, 
since  A2iP, .  -  II 21  in  this  region.  Of  course,  the 
emitted  spectrum  of  optically  thick  lines  tends  to 
be  dominated  by  surface  emission.  Nevertheless, 
one  still  might  expect  LTE  to  prevail  in  the  optical¬ 


In  Eq.  i4)  H'j2  is  the  collisional  rate  coupling 
states  i  and  i,  Pf  is  the  escape  probability  for 
the  photon  emitted  during  a  radiative  transition 
from  state  /  to  i,  is  the  spontaneous  transition 
probability,  and  U^61  is  the  radiative  recombina¬ 
tion  rate  from  the  bare  nucleus  to  NeX  n  =  l.  Equa¬ 
tion  (4)  contains  all  of  the  important  processes 
populating  and  depopulating  Nex  »  =  1  at  this  tem¬ 
perature  and  density.  The  processes  which  are 
omitted  —  collisional  recombination  lrom  NeXl  and 
its  inverse,  and  ionization  and  recombination  pro¬ 
cesses  from  and  to  Neix  —  are  relatively  unimpor¬ 
tant  compared  to  what  is  included,  and  thus  tin  ir 
omission  does  nut  quantitatively  afteet  this  analy¬ 
sis.  Equation  (4)  is  easily  solved  lor  \2  .V.  When 
the  collisional  rates  and  Einstein  ,-t's  are  inserted, 
the  ratio  ,V2  ,\;  is  found  to  be  a  function  of  the  es¬ 
cape  probabilities  for  a  particular  plasma  tem¬ 
perature.  density,  and  size.  For  Ly,  .  the  numer¬ 
ical  results  yield  ^  0.7.  and  in r  the  «=  1  re¬ 
combination  continuum  F.,.,  0.3  The  Ly.  radia- 
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Hun  leak  in  concert  with  Lyo,  y .  and  '  would  re¬ 
duce  A,  .V.  to  0.65  times  the  LTE  ratio;  inclusion 
l t  the  tree-bound  continuum  leak  brings  the  ratio 
down  to  the  actual  0.05,  which  is  0.38  that  ot  LTE. 
Photon  escape  in  Ly>  and  Lvi  makes  little  differ¬ 
ence  since  their  spontaneous  transition  probabili¬ 
ties  are  factors  of  4.4  and  13.5,  respectively,  be¬ 
low  that  of  Lyp¬ 
in  summary,  escape  of  radiation  in  the  optically 
grey  Lyp  and  NeXl  continuum  transitions  adds  to 
the  radiation  escape  in  the  line  wings  of  Lyo  and 
prevents  LTE  from  being  established  even  in  the 
plasma  interior.  While  the  strong  collisional  coup¬ 
ling  of  the  states  above  NeX  n  =  l  forces  n-  3  and 
above  into  LTE  with  n  =2,  it  in  turn  adds  to  the 
escape  of  radiation  emitted  during  decays  of  these 
levels  and  further  reduces  all  the  excited-state 
populations  relative  to  «  =  1.  Fundamentally,  the 
lack  of  detailed  balancing  of  all  of  the  NeX  radia¬ 
tion —  i.e.,  by  reabsorption— prevents  LTE  condi¬ 
tions  from  being  established. 

The  need  for  a  reduced  density  in  order  to  ex¬ 
plain  the  observed  spectrum  (from  7.8  xlO21  to 
4.5  xlO21  cm"3)  is  implied  by  two  features  of  the 
predicted  spectrum  at  7.8  xlO21  cm"3;  namely, 
the  Lyo.  Lyp,  and  Ly>  lines  are  too  wide,  and 
the  Lvo  line  exhibits  a  self- reversal  which  is  not 
observed.  We  will  first  consider  the  question  of 
the  Lyp  and  y  lines.  As  previously  mentioned,  the 
total  concentration  of  hydrogenlike  NeX  in  our 
300- eV  model  plasma  is  roughly  equal  to  the  LTE 
value.  However,  since  our  excited-state  densities 
are  much  less  than  the  LTE  amounts,  the  ground 
state  has  a  greater  than  LTE  density  —  varying 
from  75rc  greater  at  the  center  (r  =  0)  to  2.2  times 
greater  at  the  plasma  surface.  The  reason  for  this 
spatial  gradient  is  that  the  radiation  field  is  less 
intense  at  the  surface,  which  therefore,  with  less 
optical  pumping,  contains  fewer  bare  nuclei  and 
more  hydrogenlike  neon.  The  optical  depths  of 
Lyp  and  Lvy  in  collisional-radiative  equilibrium 
are  therefore  about  twice  the  previously  assumed 
LTE  values.  We  find  a  peak  optical  depth  of  2.2 
for  Lyp  in  contrast  to  the  value  of  1.2  that  was 
previously  used1  to  calculate  the  first-order  opa¬ 
city-broadening  correction.  This  results  in  a 
50  '  increase  in  linewidth  because  of  the  extra 
opacity  broadening.  Also,  w’e  do  not  find  that  the 
Ly;  line  is  optically  thin;  instead,  it  has  a  peak 
optical  depth  of  1.1.  In  summary,  significantly 
greater  opacity  broadening  for  both  Ly3  and  Ly> 
occurs  in  ORE  than  in  LTE,  which  forces  a  reduc¬ 
tion  in  density  in  order  to  reduce  the  Stark  and 
opacity  widths  to  fit  the  experimental  spectrum. 

As  noted  above,  the  300-eV  model  plasma  pre¬ 
dicts  ar.  asymmetric  self- reversal  for  Lyo.  This 
asymmetry  is  due  to  the  effects  of  plasma  motion 


The  self-reversal  itself  is  a  well-known  effect  of 
large  line  opacity  in  a  non-LTE  plasma.  Reducing 
the  linewidth  by  reducing  the  density  from  7.8  xlO21 
to  4.5  xlO21  cm"3  turns  the  self- reversal  Into  a 
"shoulder"  which  is  more  consistent  with  experi¬ 
ment.  The  increased  temperature  of  385  eV  re¬ 
duces  the  fraction  of  ground  state  NeX,  and  thus 
the  Lyo  opacity,  contributing  to  this  effect.  The 
width  of  the  opacity- broadened  Lyo  line  is  also  in 
very  good  agreement  with  experiment  when  this 
lower  density  is  assumed,  as  is  apparent  in  Fig.  2. 

IV.  EFFECTS  OF  SIMPLIFYING  ASSUMPTIONS 
ON  THE  CALCULATED  SPECTRA 

We  have  seen  above  that  neither  Voigt  profiles 
nor  LTE  are  valid  assumptions  for  the  physical 
state  of  these  high-density  neon  plasmas.  How¬ 
ever,  it  is  of  interest  to  invert  the  analysis  and 
examine  spectra  calculated  using  these  erroneous 
assumptions  order  to  estimate  the  magnitude 
of  error.  Figure  3  displays  three  spectra  which 
represent  the  results  of  employing  commonly 
used,  but  inaccurate  approximations.  The  top 
spectrum  was  computed  using  conditions  identical 
to  those  of  the  “best  fit"  to  the  experimental  spec¬ 
trum  at  the  top  of  Fig.  2,  except  that  Voigt  rather 
than  Stark  line  profiles  were  used.  Consequently, 
the  much  narrower  lines  of  this  spectrum  bear 
little  resemblance  to  the  observed  spectrum  or  to 
the  theoretical  Stark-profile  spectrum  in  the  mid¬ 
dle  of  Fig.  2.  The  bottom  two  spectra  of  Fig.  3 
were  computed  with  Stark  profiles  but  with  the 
LTE  assumption.  Since  the  true  excitation  state  of 
the  plasma  lags  behind  LTE,  when  LTE  is  as¬ 
sumed,  the  theoretical  excitation  state  is  in¬ 
creased,  resulting  in  more  ionization  and  less 
Ne  IX.  The  reduced  hehumlike  neon  is  reflected 
in  the  very  weak  NeIX  Is2-  ls3/>  ‘P  line,  which  for 
T  =  330  or  385  eV  is  much  less  intense  than  that  . 
which  is  measured.  Also,  the  NeDC  ground-stale 
recombination  continuum  is  much  reduced,  espec¬ 
ially  in  the  case  of  the  300-eV  spectrum.  Neither 
LTE  spectrum  agrees  with  experiment.  It  is  con¬ 
ceivable  that  an  LTE  spectrum  of  neon  at  a  much 
lower  temperature  might  show  a  relatively  in¬ 
creased  NeIX  lsJ-ls3p*P  line  and  thus  agree  more 
closely  with  the  experimental  spectrum  of  Fig  2. 
However,  the  lower  temperature  would  reduce  the 
collisional  rates  even  further — meaning  that  LTE 
would  be  an  even  more  inappropriate  model. 

V  SIMM  ARY  AND  CONCLUSIONS 

We  have  utilized  a  first-principles,  multistate, 
multilevel  model  of  highly  ionized,  Ingn-density 
neon  to  predict  the  emitted  spectrum  as  a  function 
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of  the  compressed  state  of  the  plasma.  The  ap¬ 
pearance  of  the  spectrum  has  been  found  to  sen¬ 
sitively  reflect  this  state,  and  the  particular  fea¬ 
tures  which  change  in  accordance  with  the  various 
plasma  properties  have  been  noted  and  the  under¬ 
lying  physics  discussed.  The  model  has  been  found 
capable  of  generating  an  excellent  detailed  fit  to 
a  previously  published  experimental  spectrum, 
and  the  fitting  process  has  led  to  a  considerable 
refinement  of  the  inferred  plasma  conditions,  in¬ 
cluding  an  estimate  of  the  extreme  sensitivity  of 
the  predicted  spectrum  to  these  conditions. 

One  should  note  that  the  authors  of  Ref.  1  at¬ 
tempted  only  to  obtain  an  approximate  temperature 
to  aid  in  inferring  the  density.  On  the  other  hand, 


it  is  our  purpose  here  to  illustrate  the  potential 
ol  non-LTE  spectroscopic  analysis  lor  obtaining 
accurate  temperatures  as  well  as  accurate  densi¬ 
ties.  We  have  shown  that  despite  the  complexity 
of  the  dynamics  of  compression,  detailed,  realis¬ 
tic,  and  accurate  modeling  of  the  radiative  output 
is  both  possible  and  fruitful. 
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ionization  from  these  levels  to  the  hydrogen-like  species.  Finally,  since  the  pumping 
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PLASMA  CONDITIONS  REQUIRED  FOR  ATTAINMENT  OF  MAXIMUM  GAIN  IN 
RESONANTLY  PHOTO-PUMPED  A1  XII  AND  Ne  IX 


I.  Introduction 

It  has  been  suggested1-4  that  population  inversions  in  plasmas  may  be 
efficiently  pumped  by  opacity  broadened  lines  from  different  ionization 
stages  of  the  same  element  or  from  different  elements  in  a  two  component 
plasma.  Experimental  evidence5’ 6  has  been  presented  for  inversions  of  the 
n  =  if-  and  n  —  3  levels  in  Mg  XII  and  Mg  XI,  which  were  pumped  by  resonant  LyO: 

and  ls2-ls2p1P  radiation  in  C  VI  and  C  V.  The  abovementioned  lasing  trans- 

c 

itions  in  Mg  lie  at  ~130  and  156  A  for  Mg  XII  and  Mg  XI,  respectively.  In 
this  paper  we  present  a  detailed  analysis  of  the  plasma  conditions  which 
would  be  needed  to  optimally  implement  two  promising  lasing  schemes  utiliz¬ 
ing  resonant  photoexcitation  with  considerably  shorter  lasing  wavelengths 

3  c 

(82  A  and  Llj-  A).  The  radiation  field-- critical  in  a  photoexcitation  process- 
is  modeled  in  detail.  The  pumped  and  pumping  transitions--as  well  as  other 
key  optically  thick  lines — are  calculated  on  a  frequency  grid  allowing  for 
accurate  modeling  of  broadening  processes  and  frequency-dependent  absorption. 
Previous  efforts  at  modeling  short -wavelength  resonantly  photo-excited  lasing 
processes  have  employed  assumed  linewidths1,  line  profiles  arising  from 
uniform  source  functions2,  assumed  power  densities3’5,  or  probability-of- 
escape  approximations4.  In  another  study',  the  pumped  plasma  was  modeled 
with  a  fine  frequency  grid  but  the  pumping  spectrum  was  assumed  to  be  a 
filtered  Planckian.  In  section  II  our  atomic  model  is  described  along 
with  the  methodology  for  its  employment  for  calculations  for  both  the  pumped 
and  pumping  plasmas.  In  section  III  the  equilibrium  results  are  presented 
for  optimum  plasma  densities,  temperatures,  and  relative  velocities  and  the 
relevant  physical  processes  controlling  these  effects  are  analyzed. 

Finally,  we  summarize  the  work  and  present  our  basic  conclusions  in  section 
IV. 


Manuscript  submitted  July  28, 1981. 
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II.  Description  of  Model  and  Calculations 


A.  Sasic  Details  of  Models 

Tne  photon  pumping  schemes  to  be  analyzed  are  the  following. 

■>  i  0 

Scheme  1  would  employ  the  Si  XIII  Is  -ls2p  P  resonance  line  at  6.650  A  to 

pump  the  A1  XII  ls^-ls3p^P  resonance  line  at  6.635  X,  creating  an  inversion 

of  the  n  *  3  and  n  *  2  levels  of  A1  XII.  Scheme  2  utilizes  the  Na  X 

ls^-ls2p^P  line  at  11.00  £  to  pump  the  Ne  IX  ls^-ls4p^P  line  at  11.001  £, 

creating  an  inversion  primarily  in  the  n  ■  3  and  n  *  2  levels  of  Ne  IX. 

Scheme  1  results  in  stimulated  emission  at  44  £  and  scheme  2  at  82  £.  Ionic 

state  and  level  densities  as  well  as  the  radiation  field  are  computed  for 

Al  using  the  model  described  in  ref.  7.  For  Si,  a  precisely  analogous  model 

to  Al — the  same  level  structure  and  transitions — is  employed.  The  Ne  atomic 

model  is  described  in  ref.  8.  This  model  possesses  an  extra  degree  of 

9 

sophistication  in  that  self-consistent  Stark  profiles  are  used  for  the  line 
opacity  rather  than  the  Voigt  profiles  employed  for  Al  and  Si.  For  Na,  no 
atonic  model  is  employed.  Rather,  the  profile  of  the  pumping  line  at  11.00  £ 
is  utilized  as  it  was  experimentally  measured  from  glass  impurities  in  a 
laser  implosion  experiment  at  the  University  of  Rochester1  .  The  multifre¬ 
quency  measured  profile  is  modified  within  the  pumped  neon  plasma  by  absorp¬ 
tion  and  re-emission  in  the  Ne  line  and  this  phenomenology  is  computed  in 
detail  using  the  flux  profile  of  the  Na  line  as  an  input  condition  on  the  Ne 
plasma.  Further  details  are  given  below  in  subsection  C.  Results  given 
below  are  calculated  for  collisional-radiative  equilibrium  (CRE). 

B.  Pumped  Plasma  Calculation 

As  previously  pointed  out^,  radiation  trapping  effects  in  the  pumped 
plasma  can  generally  be  minimized  or  eliminated  by  making  the  transverse 
dimensions  of  the  plasma  small.  We  have  followed  this  approach  in  modeling 
the  pumped  plasma — the  radiation  field  is  calculated  in  a  planar  plasma  of 
infinite  area  with  thickness  small  enough  to  insure  an  optically  thin  regime 
over  a  broad  angular  range  of  specific  intensities.  Our  objectives — given  an 
optically  thin  lasing  medium — are  to  determine  a  range  of  temperatures,  den¬ 
sities,  and,  for  scheme  1,  streaming  velocities  that  will  optimize  gain  and 
to  determine  some  of  the  tradeoffs  involved.  We  have  previously  discussed 
trapping  effects  in  the  pumped  plasma  in  refs.  4  and  7. 
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The  pumped  plasma  is  assumed  to  be  bathed  symmetrically  in  the  pumping 
radiation  which  is  calculated  (for  Si)  or  measured  (for  Na) .  The  penetration 
of  the  pumping  radiation  into  the  pumped  plasma  is  calculated  in  a  straight¬ 
forward  manner.  At  the  outer  boundaries  of  the  pumped  plasma  the  inward 
specific  intensity  along  each  ray  at  each  frequency  I  is  taken  to  be  that 
emitted  from  the  pumping  plasma  which  is  also  assumed  to  be  planar.  The 
radiative  transfer  equation  is  then  solved  in  the  pumped  plasma  with  this 

particular  boundary  condition  for  I  applied  at  each  of  the  chosen  rays  at  the 

v 

outer  boundaries  of  the  pumped  plasma.  The  calculation  in  the  pumped  plasma 
then  proceeds  by  iteration  until  steady  state  conditions  are  obtained.  Since 
fully  coupled  radiative  transport  and  rate  equations  are  solved  in  this  model, 
the  steady  state  obtained  is  a  self-consistent  collisional  radiative  equilib¬ 
rium.  The  quantity  of  primary  interest  here--the  gain  coefficient  of  the 
lasing  transition--is  obtained  as  a  linear  function  of  the  computed  densities 
of  the  upper  and  lower  states.  Finally,  the  temperature,  density  and  (for  Al) 
velocity  of  the  pumped  plasma  was  varied  in  order  to  obtain  the  functional 
dependence  of  the  gain  coefficient  on  these  quantities.  We  make  no  attempt  to 
calculate  the  depletion  of  excited  states  by  lasing,  and  thus  are  computing 
only  the  linear  amplifier  behavior. 

C.  Treatment  of  the  Pumping  Plasma 

Our  principal  objective  is  to  determine  the  effect  of  varying  condi¬ 
tions  in  the  pumped  plasma  on  the  achieved  gain.  Therefore,  we  selected  only 
one  set  of  pumping  conditions  for  each  of  the  Si  and  Na  plasmas.  Since  the 
ls2-ls2p1P  Si  XIII  and  1s2-1s3p1P  Al  XII  lines  are  off  resonance  by  0.015  A, 

the  profile  of  the  ls2-ls2p1P  line  must  be  wide  enough  to  produce  significant 

c 

intensity  at  0.015  A  (i.e.  13  Doppler  widths  at  it-OOeV)  from  line  center.  Thus 
for  a  stationary  plasma  the  line  must  be  opacity  broadened2.  If  the  resonance 
line  is  very  thick  at  line  center,  the  Lorentz  wings  will  still  exceed  optical 
depth  unity  many  Doppler  widths  from  line  center,  guaranteeing  a  wide  profile. 
Such  conditions  can  be. obtained  in  a  moderate  energy  Si  plasma  of  1.5  mm  width 
temperature  k-OOeV,  and  ion  density  8X10lscm  3.  Plasmas  similar  to  this  have 
been  realized  in  the  labor atory11.  Most  importantly,  the  CRE  calculation  for 
this  plasma  indicates  that,  depending  on  position,  62-82%  of  the  ions  are  in 
the  active  (helium-like)  stage.  Furthermore,  the  optical  depth  of  the  pumping 
resonance  line  ls2-ls2p1P  is  ~50C,  which  produces  a  very  wide  profile,  as 
shown  in  Fig.  1.  In  addition  to  the  calculated  emission  profile  the 
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intrinsic  absorption  profile  for  the  ls2-ls3piP  A1  XII  line  is  shown  for 
typical  conditions  in  the  aluminum  plasma.  In  this  case,  the  Si  resonance 
line  is  sufficiently  broadened  by  opacity  to  overcome  the  resonance  defect. 
Therefore,  under  the  reasonable  assumption  that  these  equilibrium  Si  plasma 
conditions  be  achieved  the  Si/Al  lasing  scheme  will  be  viable  at  least  in 
this  sense  that  the  resonance  defect  can  be  overcome. 

To  obtain  a  radiation  source  to  pump  Ne,  we  have  analyzed  the  spectrum 
of  a  laser- imploded  neon  filled  glass  microballoon  obtained  at  Rochester12. 
One  of  the  strongest  lines  appearing  in  this  spectrum  is  that  of  the  Na  X 
ls2-ls2p1P  line--which  arose  from  sodium  impurities  in  the  glass.  Since  we 
have  been  able  to  reproduce  the  observed  spectrum  with  a  first-principles 
non-LTE  calculation  of  the  line  and  continuum  intensities3,  the  theoretical 
calculation  which  matches  the  observed  spectrum  also  yields  the  absolute 
intensities  of  the  lines--(flux  in  ergs/cm2-sec-Hz)  at  the  outer  surface  of 
the  pellet.  Knowing  the  absolute  intensities  of  the  Ne  lines,  one  may  infer 
the  Na  line  intensity  profile  from  its  measured  intensity  relative  to  the  Ne 
lines.  Its  value  at  the  central  peak  is  6.6X103  ergs/fcm2-sec-Hz) .  It  is  this 
intensity  profile  which  we  use  as  our  pump  source  in  the  Na/Ne  calculations 
discussed  in  the  next  section. 

The  frequency  profile  is  the  one  appearing  in  the  published  spectrum-2, 
where  experimental  sources  of  broadening  are  relatively  small.  In  any  event 
deconvoluting  any  experimental  broadening  would  result  in  a  sharper  central 
peak,  which,  since  the  lines  are  perfectly  resonant,  would  give  a  better 
pump  source.  In  the  next  section  it  is  shown  that  this  experimentally 
observed  Na  X  ls2-ls2pJ'P  line  intensity  profile  is  sufficient  to  pump  a 
substantial  inversion  in  the  Ne  X  n  =  2  and  n  =  3  levels.  This  finding  is 
significant  in  light  of  the  fact  that  this  radiation  was  merely  a  consequence 
of  sodium  impurities  in  the  glass,  i.e.  no  effort  was  made  to  increase  its 
intensity  in  the  experiment.  We  note  that  only  the  radiation  of  the  Na 
resonance  line  was  used  to  pump  Ne,  whereas  the  pumping  effects  of  all  the 
calculated  Si  radiation-pumping  plus  other  lines  plus  continuum,  were 
included  in  the  A1  calculation. 


III.  Results  of  Gain  Calculations 
A.  Density  Dependence 

A  principal  question  related  to  the  time  varying  conditions  in 
the  pumped  plasma  is:  at  what  density  is  gain  maximized?  As  has  been  often 
pointed  out3’4*6’12,  there  exists  for  each  lasing  scheme  contemplated  a 
density  above  which  no  inversion  is  possible,  due  to  the  tendency  of  colli- 
sional  processes  to  bring  the  state  densities  into  LIE.  At  much  lower 
densities  one  also  expects  that  the  gain  coefficient  will  be  reduced  purely 
because  there  are  fewer  lasing  ions  in  a  given  linear  distance.  It  is 
evident  that  some  density  must  exist  at  which  the  gain  in  steady  state 
will  be  a  maximum.  This  behavior  is  shown  in  Figs.  2  and  J,  where  the 
results  of  detailed  calculations  for  various  pumped  plasma  densities  are 
presented  for  Si/Al  and  Na/Ne,  respectively.  In  each  case  the  gain  for  the 
strongest  of  the  2-3  lines  (the  2p-3d)  is  plotted  against  the  pumped 
plasma's  density.  For  Ne,  the  assumption  of  statistical  equilibrium  for  the 
n  =  3  singlet  sublevels  is  enforced,  whereas,  for  the  Si/Al  calculation, 
results  are  presented  with  this  assumption  both  enforced  and  relaxed.  The 
chosen  temperatures  approximately  correspond  to  maximum  values  for  gain, 
this  point  is  discussed  further  in  the  next  section.  Note  that  substantial 
peak  gains  -  of  1C2  cm  -1  or  more  -  are  obtained  for  both  schemes.  For  the 
Si/Al  calculations,  zero  relative  velocity  between  the  two  components  was 
assumed.  The  assumption  of  statistical  equilibrium  leads  to  an  overestimate 
of  gain  at  low  densities  ;  however,  for  densities  near  the  predicted  peak 
gain,  around  UX1C2°  cm  3,  the  overestimate  is  very  slight.  This  effect  is 
due  to  the  increasing  validity  of  the  statistical  equilibrium  assumption  at 
higher  densities  as  collisional  processes  dominate  the  populating  of  the 
sublevels. 

For  Al,  peak  gain  occurs  at  an  ion  (electron)  density  of  kXiO22 
(-i.2XlC21)  cm  3,  for  Ne  the  corresponding  figures  are  1C2C  (S.1X1020)  cm"3. 
For  hydrogenic  lasing  schemes,  Bhagavatula3’6  has  presented  reduced  variable 
equations  which  demonstrate  that  the  dependence  cf  electron  density  at  peak 
gain  on  Z  is  Z7.  We  note  with  interest  that  the  Z-dependence  of  electron 
density  at  peak  gain  implied  by  the  above  numbers  for  our  helium-like 
schemes  is  fairly  similar,  Z6'3. 
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ION  DENSITY  (10iacm"J) 

Fig.  2  —  Line  center  gain  coefficient  in  the  A1  XII  2p1P-3d1D  line  is  plotted 
vs.  A1  plasma  ion  density,  with  an  assumed  pumped  plasma  temperature  of 
100  eV.  Dual  results  for  the  assumption  of  collisional  equilibrium  (solid  line) 
between  the  3p*3d  states  and  for  a  general,  rate-by-rate  treatment  of  these  states 
are  displayed.  Characteristics  of  the  pumping  Si  plasma  are  discussed  in  the  text. 
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ION  DENSITY  (cm-3) 

Fig.  3  —  Line  center  gain  coefficient  in  the  Ne  IX  2p1P-3d1D  line  is  plotted  vs.  Ne 
plasma  ion  density,  for  an  assumed  pumped  plasma  temperature  of  65  eV. 
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B.  Temperature  Dependence 

In  Figs,  i)-  and  5  results  for  gain  vs.  temperature  are  plotted  for 
the  Si/Al  and  Na/Ne  systems  for  fixed  ion  densities  of  5XIC*9  cm  3  and 
102=  cm  3,  respectively.  The  gain  dependence  on  temperature  is  similar  to 
that  on  density  in  that  a  maximum  is  exhibited  with  a  sharp  falloff  on  one 
side  caused  by  the  variation  of  the  active  Ion  species'  populations  with 
temperature.  Note  however,  that  the  temperature  of  maximum  gain  is  much 
lower  than  one  would  expect  from  coronal  model  calculations  of  the  ionic 
species  abundances.  Such  calculations13  reveal  that  helium-like  ion  concentra 
tions  peak  at  !2CeV  and  330eV  for  Ne  and  Al,  respectively.  The  present 
calculations  predict  corresponding  peak  gains  at  ?QeV  and  lOOeV.  These 
lower  pumped  plasma  temperatures  are  necessitated  by  the  radiative  pumping 
to  the  n  =  ^  and  n  =  3  bound  levels  in  the  two  schemes  which  greatly  facil¬ 
itates  collisional  ionization.  Hence  peak  lasing  ion  abundance  is  forced 
to  much  lower  temperatures  where  fewer  electrons  are  capable  of  ionization 
from  these  bound  states.  At  still  lower  temperatures  a  sharp  gain  and 
abundance  falloff  occurs  as  the  plasma  assumes  a  more  normal  configuration 
when  the  "extra"  ionization  becomes  small.  In  these  two  lasing  systems, 
the  strength  of  the  pumping  decisively  affects  the  temperature  at  which  the 
pumped  plasma  must  be  prepared  for  maximum  gain  to  occur.  Such  effects  have 
been  noted  elsewhere6’7  In  somewhat  different  contexts. 

C.  Velocity  Dependence 

For  the  Na/Ne  system,  the  pumped  and  pumping  lines  are  within 
resonance  to  1  part  in  105,  and  thus  there  is  no  question  as  to  the  adequacy 
of  the  wavelength  coincidence.  For  Si/Al,  however,  the  wavelength  differ¬ 
ence  of  C.C15  ^  amounts  to  13  Doppler  widths  at  tOOeV.  For  a  photon  travel¬ 
ling  at  normal  incidence  to  the  pumped  plasma,  this  resonance  defect  could 
be  made  up  if  the  two  plasmas  stream  toward  each  other  at  6.8X107  cm  sec  -. 
But  given  the  wide  opacity  broadened  pumping  line  profile  (Fig.  1)  the 
functional  dependence  of  gain  on  plasma  streaming  velocity  must  be  calcula¬ 
ted  ;  this  result  is  presented  in  Fig.  6.  In  these  calculations,  the 
frequencies  of  the  radiation  incident  on  the  pumped  plasma  were  shifted 
angle-by-angle  to  reflect  the  streaming  velocities  indicated.  Even  though 
there  is  a  substantial  self-reversal  at  the  center  of  the  pumping  line, 
peak  gain  does  indeed  occur  for  the  matched  streaming  velocity  of  6.8X1 C' 
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Fig.  5  —  Line  center  gain  coefficient  in  the  Ne  IX  2p1P-3d1D  line  is  plotted 
vs.  Ne  plasma  temperature  for  a  fixed  Ne  ion  density  of  1020  cm"3. 
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Fig.  6  —  Line  center  gain  coefficient  in  the  Al  XII  2p1P-3d1D  line  is  plotted 
vs.  velocity  of  approach  of  the  Si  and  Al  plasmas.  The  Al  ion  density  is 
4  x  1020  cm*3  and  its  temperature  (the  same  as  that  of  the  pumping  Si  plasma) 
is  assumed  to  be  400  eV. 
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cm  sec-\  This  is  due  to  the  fact  that  most  of  the  pumping  radiation  is  not 
normally  incident  on  the  pumped  plasma  and  therefore  a  range  of  velocity 
shifts  are  sampled  (due  to  the  cos  3  effect)  at  any  one  physical  streaming 
velocity.  At  the  perfectly  matched  streaming  velocity  of  6.8X10^,  gain  is 
-3  times  that  of  zero  streaming  velocity  because  the  very  highest  pumping 
line  intensities  just  outside  the  self  reversed  core  are  sampled  to  the 
greatest  degree.  Having  the  two  components.  Si  and  Al,  approach  each  other 
is  therefore  helpful,  but  not  essential  to  the  scheme's  basic  viability. 

IV.  Further  Remarks  and  Conclusions 

We  have  determined  through  a  series  of  detailed  calculations  the  condi¬ 
tions  under  which  significant  gain  at  x-ray  wavelengths,  employing  the  Na/Ne 
and  Si/Al  plasmas  for  resonant  photon  pumping,  should  be  attainable  in  the 
laboratory.  Substantial  gain  at  82  X  and  44  X  for  Na/Ne  and  Si/Al,  respective¬ 
ly,  is  in  principle  achievable,  as  documented  in  Figs.  2-6.  However,  the  task 
of  setting  up  the  correct  plasma  conditions  is  not  trivial  for  a  number  of 
reasons.  For  optimum  employment  of  both  schemes,  the  temperature  of  the 
pumped  plasma  should  be  maintained  well  below  that  of  the  pumping  plasma  to 
avoid  excessive  ionization  in  the  lasing  medium.  This  could  perhaps  be 
accomplished  by  keeping  the  two  components  as  physically  separate  as  possible 
to  reduce  conductive  temperature  equilibration.  Similarly,  the  pumping  plasma 
might  be  heated  first,  and  then  the  pumped  medium  activated  through  use  of  a 
delayed  heating  pulse  or  laser  beam  to  assure  that  the  pumped  plasma  passes 

through  the  optimal  temperature  range  while  being  exposed  to  the  intense  pump- 
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ing  radiation.  Also,  in  previously  successful  experiments  ,  stepped  targets 
using  metal  plates  as  heat  sinks  have  allowed  experimenters  to  tune  the  plasma 
temperature  downward  at  certain  distances  from  the  initial  plasma  formation 
surface.  Perhaps  similar  techniques  could  be  employed  for  the  present  schemes. 
Even  though  a  lower  pumped  plasma  temperature  is  essential  for  optimum  steady 
state  gain,  substantial  gain  is  still  achievable  for  equal  pumped  and  pumping 
plasma  temperatures  (Fig.  6). 

In  the  case  of  the  Na/Ne  system,  pumping  radiation  was  generated  in  an 
actual  pellet  implosion  experiment  at  Rochester  for  which  the  ion  density  of 
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the  pellet  has  been  diagnosed  as  4.5X10  cm--3  ,  which  is  more  than  an  order  of 

magnitude  greater  than  the  neon  ion  density  required  for  maximum  gain.  In  short, 
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a  vary  dense  sodium  plasma  is  desirable  co  obtain  high  pumping  power,  but  a 
relatively  tenuous  neon  medium  is  needed  to  prevent  collisional  processes 
from  neutralizing  the  pumped  inversion.  Therefore,  a  configuration  which  is 
the  reverse  of  a  normal  pellet  suggests  itself.  One  might  compress  a  cylin¬ 
drical  glass  rod  (with  a  cylindrically  focussed  laser,  perhaps)  which  has  been 
heavily  doped  with  sodium  impurities.  This  rod  would  initially  be  encased 
in  neon,  which  would  form  a  more  tenuous  blowoff  plasma.  Or,  alternatively 
two  physically  separate  Na  and  Ne  plasmas  could  be  created  with  intensities 
and  pulse  widths  tailored  to  produce  optimum  gain  characteristics.  This 
would  certainly  allow  different  densities  to  be  produced  in  the  separate 
components ,  and  would  minimize  or  eliminate  conductive  temperature  equili¬ 
bration.  Since  ionization  of  the  pumped  plasma  by  pumping  radiation  other 
than  the  resonance  line  did  not  present  any  serious  difficulties  in  the 
Si/Al  calculation,  it  would  not  be  expected  to  present  a  problem  in  the 
Na/Ne  case.  However,  this  cannot  be  stated  with  complete  confidence  because, 
as  mentioned  above,  only  the  Na  resonance  line  radiation  is  assumed  incident 
on  the  Ne  in  our  calculations. 


ACKNOWLEDGMENT 

This  work  was  supported  in  part  by  the  Defense  Nuclear  Agency. 


C-  2  2 


REFERENCES 


1.  A.V.  Vinogradov,  I. I.  Sobelman,  and  E.A.  Yukov,  Kvant.  Electron.  (Moscow) 
2,  105(1975)  rsov.  J.  Quantum  Electron.  2,  59  (1975)]. 

2.  B.A.  Norton  and  N.J,  Peacock,  J.  Phys.  B  8,  989  (1975). 

3.  V.A.  Bhagavatula,  J.  Appl.  Phys.  L7,  L535  (1976). 

4.  J.P.  Apruzese,  J.  Davis,  and  K.G.  Whitney,  J.  Phys.  B  11.  LokJ  (1978). 

5.  V.A.  Bhagavatula,  Appl.  Phys.  Lett.  22.,  7 26  (1978). 

6.  V.A.  Bhagavatula,  IEEE  J.  Quantum  Electron.  16,  603  (I98O). 

7.  K.G.  Whitney,  J.  Davis,  and  J.P.  Apruzese,  Phys.  Rev.  A  22,  2196  (I96C). 

8.  J.P.  Apruzese,  P.C.  Kepple,  K.G.  Whitney,  J.  Davis,  and  D.  Duston, 

Phys.  Rev.  A.,  in  press. 

9.  H.R.  Griem,  M.  Blaha,  and  P.C.  Kepple,  Phys.  Rev.  A  IP,  2L21  (1979). 

10.  B.  Yaakobi,  D.  Steel,  E.  Thorsos,  A.  Hauer,  and  B.  Perry,  Phys.  Rev. 

Lett.  2°,  1526  (1977). 

11.  P.  Burkhalter,  J.  Davis,  J.  Rauch,  W.  Clark,  G.  Dahlbacka,  and  R. 
Schneider,  J.  Appl.  Phys.  ^0,  705  (1979). 

12.  K.G.  Whitney,  J.  Davis,  and  J.P.  Apruzese,  "Some  Effects  of  Radiation 
Trapping  on  Stimulated  VUV  Emission  in  Ar  XIII",  in  Cooperative  Effects 
in  Matter  and  Radiation,  edited  by  C.M.  Bowden,  D.W.  Howgate,  and  H.R. 
Robl  (Plenum,  New  York,  1977) • 

13.  V.L.  Jacobs,  J.  Davis,  J.E.  Rogerson,  and  M.  Blaha,  Astrophys.  J.  25C. 

627  (1979);  also  unpublished  calculations  for  Al. 

iL.  V.A.  Bhagavatula  and  B.  Yaakobi,  Opt.  Commun.  2k,  351  (1978). 


C  - :  5 


Appendix  D 

DESCRIPTION’  OF  WHY  RAD  CODE 

This  appendix  is  a  brief  documentation  of  the 
main  routine  and  subroutines  employed  in  the  program.  It 
is  expected  that  a  user  closely  interested  in  the  code 
will  have  access  to  a  listing;  this  documentation  should 
facilitate  understanding  of  the  code  archi tecture .  Most 
of  the  quantities  employed  in  the  input  data,  in  subroutine 
IXITAL  or  main  program  IVHYRAD,  are  described  in  comment 
blocks  immediately  adjacent  to  their  input  location. 

Figure  44  shows  a  general  flow  chart  of  the  code. 

Subroutine  IXITAL  provides  the  initialisation  of 
the  dynamic  variables  whether  it's  a  complete  job  start 
or  restart.  For  a  t  =  0  start  the  radial  distributions  of 
temperature,  adiabatic  exponent,  number  density,  magnetic 
field,  etc.  are  imputted;  a  t  >  0  restart  is  accommodated 
by  a  call  to  tape  M7.  The  number  density  initial  condi¬ 
tions  may  be  specified  by  Gaussian  distributions  of  Hj, 
polynomial  representations  of  n^,  or  general  cell -by 
cell  read  in  of  c  f=-njmj). 

Subroutine  READ1X  incorporates  the  atomic  physics 
data  for  inclusion  into  the  ATPHYS  ionization  chemistry- 
radiation  package.  Only  argon  data  for  RF.AP1N  arc 
presently  available. 
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Subroutine  ATPHYS  i'eeds  the  ground  plus  first 


excited  state  physics  into  the  hydrodynamic  cells  at  each 
time  step:  within  this  routine  estimates  of  ionicatien 
level,  adiabatic  exponent,  and  the  magnitudes  of  various 
radiation  forms  are  made.  The  flag  KRAD  =  1  selects 
this  routine. 

Subroutine  YOI.TG  incorporates  the  voltage-time 
profile  into  the  model;  linear  interpolations  between  a 


maximum  of  50  normal  iced  time  points  (t/t 


pulse 


i  are 


employed  as  the  representations. 

Subroutine  CURG  incorporates  a  current-time  profile, 
if  desired  instead  of  voltage,  into  the  model.  Linear 
interpolations  between  a  maximum  of  50  normal  iced  t.imc 
points  are  employed.  The  flag  1CUR  =  0  chooses  the  voltage 
input,  1CUR  •=  1  chooses  the  current  input. 

Program  IVHYRAP  is  the  main  routine  in  the  code. 

It  accepts  the  namelist  REAP  I \  of  many  of  the  control 
variables,  performs  the  main  solution  operations,  specifically 
the  extenal  circuit  coupling  and  the  calling  cf  the  1  1  HI G ! 
numerical  integration  algorithms,  specifies  the  time  stop 
loop,  prints  out  the  results  and  sets  up,  if  desired,  the 
information  to  be  stored  on  file  for  the  next  restart. 

Subroutine  RTPIIYS  serves  a  similar  function  that 
ATPHYS  does;  it  translates  the  i on  i  ca t i on  -  rad i a t i on  output 
from  the  multicell  real-time  package  or  the  TLRRYS  curve-fit 
package  to  conform  as  an  input  to  Kin  RAD.  If  the  ILRRYS 
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radiation  package  is  employed  then  ionization  number  and 
total  chemical  potential  (used  for  the  calculation  of  the 
isentropic  exponent)  are  also  determined  within  RTPHYS. 

In  addition,  a  benchmark  radiationless  hydrodynamic  option 
(  I II Y D  >  0.)  may  be  called  within  this  subroutine. 

Subroutine  TERRYS  links  together  the  curve  fit 
radiation  package  and  writes  out  the  total  (line 
plus  continuum)  radiation  power  at  each  cell  (erg/sec  cniJ  .  , 
as  well  as  the  cross  sectional  area  integrated  power 
lerg/sec  cm)  for  the  total  amount,  the  line  radiation  <  1  keY 
photons  ('L-shelli,  >  1  keY  (K-shell),  and  continuum. 

Subroutine  DU ST SR  incorporates  the  curve  fits  of 
the  aluminum  or  argon  MCELL  optically  thin  radiation  pack¬ 
age,  with  polynomial  dependence  of  volume  emission  rate 

s  l  d  -  s 

(erg/sec  cm')  on  electron  temperature  at  n^  =  10  '  cm  ' ; 

dependence  at  other  densities  is  estimated  by  scaling 

estimates  of  power  law  form. 

Real  function  RADSSP  computes  the  opacities  from 
the  radiating  cell  out  to  the  plasma  boundary  along  the 
photon  path  length;  free-bound,  K  and  L  shell  line  and  a 
single  average  frequency  free- free  radiation  modes  are 
considered.  If  IBCHOC  =  1  a  more  exact  ten  frequency  IT 
Bremsstrahl ung  representation  is  employed. 

Subroutine  FI. BAR  incorporates  the  computation  of 
the  ten  frequency  FF  path  length  along  two  rays  emitting 
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from  the  cell  of  interest.  Integration  is  carried  through 
the  range  of  number  densities,  ionication  levels  and 
temperatures  along  the  path. 


Real  function  RCSSIG  evaluates  the  recombination 
continuum  absorption  cross  section  needed  for  evaluation 
of  the  bound- free  optical  thickness  in  RADSSP. 

Real  function  PESFBC  computes  the  probability  of 
escape  for  bound-free  radiation,  while  real  function 
PESLKT  computes  the  probability  of  escape  for  line 
radiation. 

Real  function  LBAR  computes  the  path  lengths  for 
line  and  bound- free  radiation  transport,  with  integration 
carried  through  the  range  of  number  densities  and 
temperatures  along  the  path. 

Subroutine  DTSET  determines  the  integration  time 
stop  in  terms  of  the  collapse  dynamics  of  the  discharge. 
The  minimum  time  step  over  the  radial  cells  is  evaluated 
on  the  basis  of  the  Courant-l-redericks-Lewy  sound  speed 
condition,  the  ratio  of  the  thermal  energy  to  the  loule 
heating,  the  ratio  of  the  thermal  energy  to  the  local 
radiative  transport,  and  the  ratio  of  the  thermal  energy 
to  the  previous  time  step  internal  energy  change.  The 
minimum  of  these  results  is  taken  and  employed.  The  last 
criterion,  based  on  the  internal  energy  change,  is  usually 


the  most  restrictive.  The  Alfven  speed  is  also  determined 
here  to  establish  minimum  inertia  and  density  in  the  outer 
regions  of  the  implosion. 

Subroutine  DRSET  evaluates  the  dynamic  rezoning 
cell  widths  based  on  the  local  cell  ion  number  densitv. 

IDR  s  1  selects  uniform  incrementation,  IDR  ^  2  imposes 
the  rezoning.  The  cell  center  radii  and  the  annular 
incremental  volume  are  also  calculated  here. 

Subroutine  BSOLVE  incorporates  the  tridiagonal 
algorithm  for  the  solution  of  the  self-induced  magnetic 
field  as  a  function  of  the  external  c i rcuit - imposed 
boundary  condition  and  the  assumed  resistivity  function. 

Subroutine  TRANS  specifies  the  transport . coef fic ient 
of  thermal  conductivity,  electron-ion  temperature  relaxa¬ 
tion,  thermoelectric  effect  and,  in  particular,  the  anomalou 
resistivity.  The  choice  of  classical  resistivity,  core¬ 
corona  anomalous  resistivity,  or  one  based  on  the  ion 
sound  speed  is  made  with  flag  1TURB  10,  !,  2,  respectively1. 

Subroutine  THIFT  involves  the  tridiagonal  algorithm 
for  the  fully  implicit  solution  of  the  energy  equation 
thermal  diffusion  source  term.  This  calculates  the  final 
temperature  value  within  a  time  step  to  represent  the 
complete  electron  temperature  and  total  internal  energy 
conservation  equations. 


Subroutine  EQU1L  evaluate?  the  effects  of  energy 
equilibration  between  the  electrons  an!  ions.  Also,  for 
comparative  purposes,  black  body  heat  flux  is  evaluated  here. 

Subroutine  ECHEK  contains  the  evaluation  of  the 
local  energy  balance  between  input  electrical  energy, 
hydrodynamic  internal  energy,  and  radiated  energy. 

Subroutine  ROSSL  contains  a  double  linear  interpola¬ 
tion  (density  and  electron  temperature!,  giving  the  Rosseland 
coefficients  for  the  black  body  radiation  calculation. 

The  ETBFCT  routines  (ETBFCT,  YELOCE,  NCR  IDE, 

SOURCD,  0GR1DE,  COXSRE J  are  exactly  as  found  in  Ref.  .1 
except  some  of  the  comment  statements  were  not  punched. 

The  MCELL  multicell  routines  specifically  called 
in  KHYRAD-MCRAT ,  UTIL  IT,  and  PRRAP  connect  the  hydro 
package  to  the  full  multicell  chemistry  plus  radiation 
transport  package.  The  calling  of  these  routines  has  to 
be  specified  in  .JSL  statements. 

The  XRLLIB  qualitative  plotting  code  is  employed 
through  the  calling  of  ASXASY  and  IPLOT. 

The  input  cards  are  first  a  XAMEL1ST  group  for 
a  WHYRAD  input,  then  another  group  of  standard  input  for 
IXITA1. .  Then  [for  the  TERRY S  package  choice,  employed 
in  the  example  calculations  of  Section  Yi  the  Bremsstrahlung 
radiation  frequency  table  is  read  in,  in  IXI1AL,  for  use 
in  subroutine  RTPHYS.  Initial  conditions  over  the  radial 


D-b 


distribution?  come  next  (  read  in  INITALt  starting  from 
RAD.N  ( I  l  through  BO  i  I  )  {magnetic  field}.  For  a  restart, 
quantities  from  EMAS,  EADP,  ...  (statement  505)  through 
to  READ  (M.105)  ( BO ( 1 i ,  1=1, \Rj,  statement  161,  are  input 
from  the  storage  file.  Finally,  the  voltage  distribution 
is  input  through  YOLTG. 
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Appendix  E 
LI  SI  Oi-  SYMBOLS 


a  B  z 

b  Coefficient  in  optically  thin  represent  at 

B  Azimuthal  magnetic  field 

B  Gaussian  density  distribution  coefficient 

c  Speed  of  light 

C^  Ion  sound  speed 

C  Sound  speed 

P  Source  function  term 

e  Electron  charge 

E  Plasma  internal  energy 

Iotal  chemical  potential  energy 
h  Planck's  constant 

!  iotal  axial  current 

IP  Ionization  potential 

j  Axial  current 

k  Boltzmann's  constant 

k  lhermal  conductivity 

Radiation  path  length 

I  Electrode  spacing  discharge  length 

!  C i rcu  i  t  i nductancc 

m  Specie  mass 


ion 


roue  i'  in  re  r  on  i  ny  algorithm 
Pinch  mass 


'! 

M 


n 


\R 


P 


P 


\1> 


r 


'  v. 


K 

Ry 

t 

T 

U. 

v 

V 


-  1 
:  \ 


Number  density 

Coefficient  in  optically  thin  representation 
Computational  cell  number 
Pressu re 

Probability  of  escape  function 

Radiative  power 

Thermal  flux 

iieatiny  term 

Radial  coordinate 

Return  circuit  radius 

Gas  constant 

Pinch  radius 

Ry dbe r  y 

1  ime 

1 empe ra t  ure 

Axial  current  velocity  (drift  speed; 

Radial  velocity 
briviny  volt  aye 
Plasma  volt aye 
Axial  coordinate 
Circuit  impedance 
1  on  i  cat  ion  number 
Atomic  number 


Ihermal  flux  coefficient 


isentropic  exponent 
Axial  electric  field 
fm i i on  coefficient 
Res i s t i vi ty 

Absorption  coefficient 
Coulomb  logarithm 
Radiation  frequency 
Dens  i  t  y 
Cross  section 
Relaxation  time 
Optical  depths 

General  iced  density  function 
Cyclotron  frequency 
Plasma  frequency 


c 

rr 

I 


) 

K 


Bound- free  radiation 
electron 

Free- free  ( Bremsst rah  lung l  radiation 
I  on 

Radiative  frequency  integer 
K- shell 
long  path 


L- shell 


s 


At  outer  pinch  radius 
Short  path 
0  At  pinch  centerline 

an  Anomalous  transport  coefficient 

ci  Classical  microscopic  transport  coeff 

O'!  Optically  Thin 


c  i  en  t 
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Appendix  1-' 

PROMINENT  FORTRAN  VARIABLES 


FORTRAN  SYMBOL 
AJ 1 1  I  ) 

AXE  t  I  ) 

AN  1(1) 

BN  (  I  .) 

BO  (  I  l 
DR  (  I  ) 

EDOTd (I  ) 
EFIELDl 1  I 
EXEVV  (  I  t 
ETA  l  I  ) 

FLPLAS 
FYOLl  I  ! 

GAMMA (  I  t 
1’NEtV  (  I  ) 

PRAD  ( 1  ) 

RADN (  I  ! 

RB I  I  l 
RC  (  1  ) 

R.NEIVI  I  1 


VARIABLE 


B  field  (new) 
B  field  (old) 
dr 
c  i 


c 

1. 

cell  volume 
Y 
P 

p 

1  RAD 

cell  center  r 
cell  boundary  radiu 
cell  center  r 


RVOLDl I ) 


